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For aggregated combat simulation models, the methods for calculating force 
attrition must be based upon sound mathematical formulations and parameter esti- 
mations. With an inherent lack of representative combat data for modern warfare 
scenarios, one effective method for determining the required parameter estimates is 
to thoroughly analyze the output from a stochastically based high-resolution combat 
model. It is this development of attrition parameters process, which so profoundly 
influences the validity of aggregated simulations, that lacks any comprehensive doc- 
umentation or mathematical justification within the modeling community. By ex- 
amining the development and validity of these processes for parameter estimation, 
valid attrition calibration formulae can be determined and used within force attri- 
tion algorithms in order to more precisely and justifiably model aggregated combat 
operations. The establishment of a user-friendly test bed for examining this attri- 
tion rate development process will play a major role in solidifying the understanding, 
implementation, and validation of current and future process techniques. 
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I. 



INTRODUCTION 



Combat modeling and simulation are important analysis tools available to mil- 
itary leaders and their supporting organizations which provide valuable insight into 
the study of military operations across the spectrum of conflict. Combat models 
are employed frequently to conduct investigative studies ranging from peacekeeping 
operations, which are relatively narrow in scope, to conventional maneuver warfare 
executed at the nation-state level. Commanders and their supporting planning staffs 
regularly seek the resulting output from the execution of combat models and simu- 
lations. Prior to the employment of man and machine into a hostile environment, 
tactical commanders have the capability to develop and assess their maneuver plans 
by use of combat models and simulations. Additionally, combat models and simu- 
lations are employed as training tools to rehearse command and control actions and 
operational procedures. While models and simulations are not in themselves predic- 
tors of the outcomes to military operations or plans, their ability to provide samples 
about possible outcomes is welcomed as a preliminary step prior to placing soldiers 
and equipment in harm’s way. As such, the designers of these models constantly 
attempt to refine the representation of warfare to ensure the model is as accurate 
and believable as possible. It is important to note that many of the interactions and 
human factors which are inherent to a military force are difficult to replicate within 
a combat model or simulation. This fact must be recognized by those employing 
combat models and simulations for mission analysis and maneuver warfare planning. 

One of the most important capabilities of combat models and simulations is 
their ability to calculate the attrition of personnel and systems of the engaged military 
forces. The prediction, analysis, and management of the attrition to personnel and 
equipment in combat are highly regarded as important tasks by military leaders and 
their staffs. And, as a result, the ability of combat models and simulations to capture 
attrition data is often a measurement of the model’s usefulness. What follows is a 
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discussion of the development of attrition measurements within the context of combat 
models and simulations. 

A. BACKGROUND 

Several important principles and concepts underlying the general class descrip- 
tion of combat models and simulations are important to this development of attrition 
measurement. Combat models and simulations are generally subdivided into two 
categories, namely, high resolution and aggregated. 

High resolution models are generally associated with smaller force compo- 
sitions where each combatant and its associated weapon or maneuver system are 
explicitly represented in the simulation. The complicated combat processes are de- 
composed into highly detailed sequences of events and activities which are carried 
out explicitly within the simulation. This detailed sequence of events and activities 
lends itself to simulation time steps which are very short. It is not uncommon for 
high resolution models to represent time explicitly to the second, updating impor- 
tant information almost continuously. The goal of the high resolution model is to 
accurately model each combat process so that the simulation results are traceable 
using mathematical formulas and constructs. The characteristics associated with the 
entities represented in the simulation are related to their engineering specifications, 
physical performance data, and specified behavioral assumptions. The existence of 
this information link from the engineering level through the behavioral engagement 
assumptions is the single most important quality of a high resolution model. The 
capability of representing weapon systems and equipment at the engineering level 
enables an analyst to observe the possible effects of alternative system specifications 
on a pre-defined simulation. Individual system contributions to the overall simulation 
are effectively analyzed in high resolution models. Another feature associated with 
high resolution models is based upon their reliance on stochastic processes. Nearly 
all of the simulation algorithms above the engineering specification level are proba- 



2 



bilistic in nature. For example, consider a firing sequence. The firing system goes 
through several sequential actions prior to firing its weapon. These include target 
acquisition, target selection, firing, impact, and damage assessment. As these events 
are not of a fixed length in actual combat, their representation in simulation relies 
upon the use of highly probabilistic random variables. Calculations of this nature 
can lead to variance among independent trials of a given combat simulation. In order 
to obtain a reliable result with respect to the variance between stochastic calcula- 
tions, the modeler must execute multiple runs of the same simulation to obtain the 
“averaged” result. One of the most significant limitations of high resolution models 
is their sheer size and complexity. The computer codes supporting high resolution 
models are highly complex and often quite large, stretching hardware and software to 
their limit. These computer codes are expensive to develop, run, and maintain. The 
complexity associated with representing each individual entity of the combat scenario 
at the engineering specification level causes another significant limitation of high res- 
olution models, namely, force size. High resolution models are limited to relatively 
small forces. In ground maneuver terms, this limitation is normally associated with 
an upper bound of roughly 1000 combatants [Ref. 1] or on the order of battalion to 
brigade levels of combat. 

Aggregated models, as the word implies, maintain varying levels of aggregation 
of the represented combatants. For example, an aggregated model may be capable 
of representing a group of riflemen with an associated consolidated firepower capa- 
bility rather than the individual soldiers. The model then portrays the aggregated 
capabilities of the group of riflemen represented graphically as a single entity. There- 
fore, an entity representing a platoon of 30 riflemen would maintain the firepower 
and maneuver capabilities associated with the 30 individual weapons systems carried 
by the 30 riflemen. The extent of aggregation within combat models varies based 
upon the level of combat being modeled and the intended use of the model. The 
aggregation process is natural within the confines of military combat models because 
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of the hierarchical command structure inherent to maneuver forces. For example, 
the natural aggregation process associated with a generic ground maneuver organiza- 
tion is individual combatant, squad, platoon, company, battalion, brigade, division, 
corps, army, and theater. Depending upon the use of the model in question, one 
of these natural groupings (normally no lower than brigade) would mark the lowest 
level of representation within the simulation. Individual combatants and groups of 
combatants which are smaller than the basic size are not explicitly represented and 
the attributes associated with these combatants are lost to the simulation. Two im- 
portant and distinctive aggregation patterns associated with models of this type are 
homogeneous and heterogeneous aggregation, and, the distinction between them is 
particularly important for attrition calculations. A homogeneous aggregation is one 
where the combat power of a unit with multiple and distinct combatants is combined 
into a single measure (or perhaps one for ground combat systems and a second for 
aircraft systems). Attrition computations are then based upon the relative power of 
the aggregated forces in battle subject to their combat power indices. Contrastingly, 
heterogeneous aggregation maintains a count of the number of surviving weapons sys- 
tems of each distinct type. This method of aggregation allows the model to represent 
the different levels of effectiveness of a particular weapon system against its possible 
targets. This form of aggregation permits more accurate attrition modeling than the 
homogeneous model resulting in less information loss. As a result, the trend of most 
modern aggregated models is to follow a heterogeneous aggregation process. 

As the size of the modeled entity increases in aggregated models, the preci- 
sion to which they can be located on the simulated battlefield may be decreased. A 
modeled combatant may maintain a center location and a frontal width where the 
individual entities making up the aggregated unit are assumed to be somewhere in the 
unit region. Information about the individual entities within the aggregated unit may 
be kept in the form of averages, but any specificity about individual locations is lost. 
Thus, the first significant consequence of aggregation is the loss of individual char- 
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acteristics within aggregated units resulting in the use of their averaged properties 
where any inherent fluctuations have been smoothed. Analogous to the smoothing 
of individual variations is the handling of stochastic processes. Since aggregate mod- 
els do not maintain a reference to what each individual combatant is doing at any 
given time, the complex probabilistic events found within high resolution models are 
replaced by averaged behaviors. The firing process described earlier (target acquisi- 
tion, target selection, firing, impact, and damage assessment) within the context of 
high resolution models may be replaced in an aggregated model by computing the 
rate at which an average tank kills enemy tanks. This derived attrition rate implic- 
itly defines each of the probabilistic processes described above with the stochastic 
variability suppressed. This fact leads to a descriptor generally associated with ag- 
gregated models, namely, deterministic. With the engineering specification level of 
detail and the probabilistic processes being suppressed within the aggregated model, 
many assumptions and scenario details get absorbed by the attrition rate number, 
which in itself is a highly variable entity. 

Lastly, aggregated models lose the information associated with event sequenc- 
ing since they do not keep track of individual actions. Thus, the precise times of 
critical events, such as target kills, are not available from aggregated model output. 
This results in the process averages being processed over relatively long periods of 
time (ranging from a few minutes to a 24 hour period) where the computation of total 
casualties occurs at the end of the time period. Again, careful consideration must be 
given to the attrition rate formulation in that there are at least two sides involved in 
the exchange of fires. Not all combatants on either side will survive to continue firing 
for the duration of the time period. For example, consider the following example and 
resulting conclusions: 

Extreme care must be taken to make the underlying assumptions clear 
in computing the rate data for such models. For example, a daily attrition rate 
for aircraft stationed at an airbase which is under attack may be inconsistent 
with the day’s plan for those aircraft (they may be flying missions and thus 
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not on the ground when their base is attacked). A simulation which assesses 
attrition on an aggregated daily cycle cannot know whether the airbase attack 
coincided with aircraft on the ground, so an assumption must be made and 
figured into the attrition rate data. 

Overall, then, the outstanding feature of aggregation in combat models 
is information loss . As a result, it is necessary to model combat processes using 
average behavior rather than individual behavior. Although values for the 
average process rates are generally computed outside the combat model, the 
complex nature of these rates demands that we devote considerable effort. ..to 
explaining their derivation. [Ref. 1] 

An obvious question with regards to the given strengths and limitations of the 
described combat modeling and simulation approaches is why high resolution models 
are not extended to larger force sizes thus maintaining the admirable qualities of the 
high resolution definition. Unfortunately, the difficulty in making this extension lies 
within the limitations of computer hardware and software and the expansive execu- 
tion time which would be associated with such large scale high resolution models. 
The sheer number and complexity of the individual combatants and their associated 
characteristics render the use of high resolution combat models for all combat simu- 
lations impractical and, in some cases, impossible. Additionally, the complex nature 
of command, control, and communications (C 3 ) between echeloned organizations is 
extremely difficult to model. The synergism and inherent lack of measurability of 
these traits would require many assumptions on the part of the modeler possibly re- 
ducing the degree of model believability. Thus, the resulting segregation of combat 
models into the high resolution and aggregated categories is a matter of necessity. 
It is important to note, however, that the two categories may not be entirely inde- 
pendent during their execution due to the aggregated models’ possible reliance upon 
the outcomes of high resolution battle simulations for the derivation of attrition rate 
parameters. While there are several currently employed methods for calculating force 
attrition within aggregated models (to be presented in Chapter III), it remains true 
that the attrition of personnel, equipment, and supplies resulting from an engagement 
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must, at some point, be calculated within an aggregated simulation. The compromise 
which is of specific interest to this research is the algorithmic application of mathe- 
matical constructs to determine aggregated force attrition based upon high resolution 
results. For example, when an aggregated combat simulation reaches a point where 
two opposing forces engage in battle, the model may implicitly rely upon high res- 
olution simulation results generated under similar environmental, force composition, 
and tactical conditions. The end result of the battle executed within the aggregated 
simulation would be new force levels for each of the engaged combatants derived pri- 
marily from a deterministic replay of a representative battle using the attrition rate 
coefficients derived from the high resolution simulation output. The intermediate 
process of estimating the attrition rates and replaying the battle mathematically is 
the conceptual link between high resolution and aggregated models. With the line di- 
viding high resolution and aggregated models clear in definition only, both categories 
have an associated collection of applications which, when combined, adequately meet 
the needs of combat modeling and simulation. 

B. APPLICABILITY TO CURRENT PROCEDURES 

In the United States Army and, moreover, the joint warfare community, one 
of the most frequently used high resolution models which provides the data necessary 
to generate the attrition rates for use in aggregated models is called the Combat 
Sample Generator (COSAGE) model. The COSAGE model is a (fairly) high reso- 
lution, stochastic model of division level ground combat which maintains a combat- 
ant resolution to friendly platoons and enemy companies. Within each maneuver 
unit, a heterogeneous list of individual weapon systems is arranged in combat forma- 
tions, interactions between opposing weapon systems are computed, and individual 
weapons system survival is governed by the laws of probability. With respect to 
aircraft, close air support of the ground battle by tactical aircraft is resolved down 
to a detailed simulation of the flight of individual aircraft. At a point in the exe- 
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cution of an aggregated model where two opposing ground forces are positioned to 
engage in combat, the representative battle simulation is carried out as a mathe- 
matical algorithm based primarily upon the results of a high resolution simulation 
(COSAGE). The killer/victim results and the associated time series of individual 
attrition are analyzed and mathematically manipulated resulting in the development 
of attrition rate estimates for each individual weapon systems of each engaged force. 
This process is referred to as the attrition calibration (ATCAL) Phase I process and 
will be discussed in detail in Chapter IV. These derived attrition rates are then 
applied to a deterministic algorithm which, when executed or “replayed” with the 
force levels present in the aggregated model, calculates resulting force levels (by 
system). These updated quantities are then passed back to the aggregated model 
as it continues through its remaining processes. This secondary “replay” process 
is commonly referred to as ATCAL Phase II. Currently, several theater level mod- 
els have the capacity to employ ATCAL Phases I and II. These include Tac tical 
Warf are (TACWAR), Concepts Evaluation Model (CEM), and Integrated Theater 
Engagement Model (ITEM). Whether a particular theater level aggregated model 
employs the ATCAL algorithm as its method for calculating attrition is largely de- 
pendent upon the respective developers’ and users’ understanding of ATCAL’s theo- 
retical development and mathematical principles, which, in themselves, remain to be 
areas of active research. See [Ref. 2]. 

C. PROBLEM STATEMENT 

In order to support current and future investigations into the theoretical un- 
derpinnings and mathematical stability of the defining algorithms within ATCAL 
Phases I and II, researchers must have a “user friendly” method of obtaining and 
analyzing the necessary high resolution output data and validating the conjectures 
present in said algorithms. Additionally, as a starting point for the mathematical 
analysis, researchers must have a thorough understanding of the basic principles as- 
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sociated with the various methods of calculating force attrition. This capability for 
analysis and foundation for understanding do not exist outside of the United States 
Army Concepts Analysis Agency (CAA) and its proprietary implementation of the 
COSAGE model and ATCAL Phase I and II algorithms. The intent of this thesis is 
to develop the framework for a test bed capable of conducting detailed analysis of the 
attrition problem. By developing realistic high resolution ground combat simulations 
using the Janus simulation software, capturing the detailed killer/victim results and 
the associated time series of casualties, and applying these results to the study of the 
aggregated attrition problem, continued and specialized future research can be con- 
ducted in this discipline. Specifically, as an immediate result of this research, future 
inquires can be made into the validity and appropriateness of the currently accepted 
versions of the ATCAL Phases I and II algorithms aiding in the proper documenta- 
tion of their conjectures and results. Moreover, with this vehicle for comparative 
analysis, efforts can expand to the development of alternative attrition algorithms to 
meet the needs of future aggregated modeling efforts. 

D. JANUS SIMULATION SOFTWARE 

As the initial step in creating a functional test bed for the study of attri- 
tion, one must choose a vehicle to conduct the sample data base of high resolution 
simulations. Janus is the primary high resolution, interactive wargaming simulation 
software used by the United States Army for the modeling of brigade sized (and be- 
low) operations, and, its widespread availability within the military community, ease 
of use, and verified credibility make it an obvious choice for meeting the demands of 
this proposed analysis. Named for the two-faced Roman god who was the guardian 
of portals and the patron of beginnings and endings, Janus had its beginnings as a 
nuclear effects modeling simulation developed by the Lawrence Livermore National 
Laboratory of the University of California. [Ref. 3] Its applicability to tactical train- 
ing and operations analysis quickly made it functional to the United States Army 
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Training and Doctrine Command (TRADOC), Training Analysis Center (TRAC). 
The current version is intended to satisfy both the combat development and training 
communities. Its versatility allows it to be used interactively as a man-in-the-loop 
simulation to replicate the effects of command and control on a realistic battlefield 
during the simulation as enemy activity develops. As an analysis tool, Janus can 
be used non-interactively thus reducing the variability introduced by human decision 
making during the simulation execution. Janus is a stochastic simulation in that the 
system determines the results of actions (like direct fire engagements) according to 
laws of probability. While the principle focus of Janus is on ground maneuver and 
artillery units, it also models weather and its effects, day and night visibility, engineer 
support, minefield employment and breaching, rotary and fixed wing aircraft, resup- 
ply operations, and chemical environments. The inherent flexibility of this simulation 
software allows its user to accurately model a multitude of combat scenarios. 

E. SCOPE AND LIMITATIONS 

This research is unique in that there has been no other attempt to establish 
an independent capability for analyzing the aggregated attrition problem, specifically 
the ATCAL Phase I and II algorithms. This work will validate the capability of con- 
structing a high resolution ground combat simulation in Janus for the purpose of 
deriving attrition rate parameters using the ATCAL Phase I and II algorithms. Ad- 
ditionally, there has been no attempt to configure the Janus output data to conform 
to the requirements of the current ATCAL Phase I and II algorithms (currently exist- 
ing in the FORTRAN computer language). This procedure and the transformation of 
the output data is practical only within the confines of analysis and algorithm valida- 
tion because of the time expenditure required to develop the ground combat scenario 
in Janus and the nonexistence of automated software to link the Janus killer/victim 
results to current aggregated models. With initial efforts confined to the manual com- 
pilation of Janus output data and the explicit formulation and calculation of basic 
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mathematical models, this research could be expanded to focus on the compatibility 
and automation of the data linkage between Janus the current and/or future ATCAL 
computer codes. However, this is not the focus of this thesis. And, as only one specific 
ground combat simulation will be presented as part of this thesis, follow-on work is 
applicable in the development of additional scenarios with varying mission objectives, 
force definitions, environmental conditions, and terrain selections in order to general- 
ize results otherwise limited to specific simulation parameters. A robust collection of 
simulations would be required to properly validate the conjectures presented in the 
problem statement. As such, the results presented in this work relate only to the 
specific scenario and situations modeled and should not be considered to be the final 
solution. 
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II. 



THE SIMULATION SCENARIO 



In order to generate the required high resolution output data with Janus for 
the development of attrition rate numbers using the ATCAL Phase I and II algo- 
rithms (analogous to the process described in Chapter I, Section B), the author has 
chosen to simulate a heavy brigade combat team (BCT) combat scenario. The typical 
heavy BCT includes a heterogeneous mix of direct fire, indirect fire, and supporting 
maneuver systems. Generically, the force would include tanks, armored personnel car- 
riers, artillery systems, air defense artillery systems, engineering equipment, wheeled 
vehicles, air assets, and obstacles. The BCT is the preferred level of combatant or- 
ganization for this study because of the described combined arms composition and 
the ability to portray this force in Janus. As a high resolution model, Janus can 
effectively support a force type of this size and capture the detailed results of all 
direct and indirect fire engagements. With these results, the attrition rate estimation 
process (ATCAL Phase I) and the deterministic battle replay (ATCAL Phase II) can 
be conducted. 

The tactical scenario chosen for this combat simulation is one of practical 
relevance. As this research focuses primarily on the modeling of ground combat 
as part of an overlying larger tactical operation, the author has chosen to simulate 
the mechanized task force operations in desert environmental conditions as a base- 
line scenario. As the Janus simulation package can be an interactive training tool 
when employed as a man-in-the-loop simulation, it is often used by TRADOC at the 
United States Army’s combat arms training centers. Specifically, the Armor Officer 
Advanced Course (AOAC) located at the Unites States Army Armor Center, Fort 
Knox, Kentucky, employs Janus as a tool to train and evaluate the mission planning 
and tactical decision making skills of its officer students. With a specified mission of 
training the current and future leaders of the United States maneuver warfare orga- 
nizations, AOAC and its focus on modern ground maneuver tactics, techniques, and 
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procedures provide an excellent simulation opportunity for this research. The Janus 
simulation used during this research was executed by AOAC Class 97-3 at the Fort 
Knox Brigade Battle Simulation Center. [Ref. 4] And, the simulation is available 
for review through coordination with Professor Bard K. Mansager, Department of 
Mathematics, Naval Postgraduate School. The stochastic nature of the simulation is 
defined not only by the probabilistic events which occur as a function of the intrinsic 
algorithms, but also by the human factors associated with the particular officer stu- 
dents who executed the simulation. If multiple runs of this scenario were executed, 
each would produce different results based upon these stochastic processes. 

A. MECHANICS OF THE SIMULATION 

At the highest level of organization, Janus allows the definition of six combat- 
ant sides where each side could represents an autonomous force. Under each side , 
forces can be further organized into task forces where modeled weapons and maneu- 
ver systems can be organized and manipulated in groupings that mirror their real life 
task organization. After the defined combatants are organized into task forces, one or 
more of these task forces are output to a specific computer workstation for user con- 
trol. Janus permits up to 24 computer workstations to be in use during a simulation. 
The 24 workstations are divided into side 1 through side 6 terminals based upon the 
number of users assigned to manipulate the respective sides and task forces. [Ref. 5] 
The simulation used for this research makes use of sides 1 and 2 and will hereafter 
be referred to as BLUE and RED, respectively. This assignment represents the roles 
of the combatants during the simulation. In this simulation, BLUE represents the 
friendly Unites States forces being controlled by the AOAC 97-3 student officers and 
RED represents a Soviet-style enemy forces being controlled by designated opposing 
force personnel. 
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1. BLUE Force Composition and Organization 

The BLUE force (BLUEFOR) is made up of the United States maneuver 
and weapons systems appearing in Appendix A, Table XIII on page 64; and, when 
organized for combat, these systems represents a heavy BCT. The BLUEFOR is 
considered to be “balanced” in terms of its armor to mechanized infantry ratio. As 
an integrated combined arms fighting force, the BLUEFOR systems are organized into 
a typical BCT hierarchical command and control structure as detailed by Table I. For 
the purposes of the AO AC student officers who executed the Janus simulation, this 
fictitious force was designated as the 1st Brigade Combat Team (1 BCT). 

Unit Description Major Systems 

4 Direct Fire Maneuver Battalions MlAl, M2 

2 Indirect Fire Artillery Battalions M109A6 



1 Air Defense Artillery Company BSFV 
1 Chemical Support Company IF TRU 
1 Engineer Support Battalion AVLM 

1 Logistics Support Battalion IF TRU 

1 Tactical Air Control Party A10 (6) 



Table I. BLUEFOR Task Organization for AOAC 97-3 Janus Simulation 
[Ref. 6] 

2. RED Force Composition and Organization 

The RED force (REDFOR) is made up of Soviet style equipment listed in 
Appendix A, Table XII on page 63; and, when organized for combat, it represents 
a less than full strength Soviet style motorized rifle regiment [MRR(-)]. Note the 
greater than three to one ratio of mechanized infantry to armor, implying that the 
REDFOR is better equipped for mechanized infantry operations than tank warfare. 
Additionally, the REDFOR has a significant number of indirect fire artillery systems 
at its disposal. The defensive REDFOR MRR(-) is organized within the hierarchical 
command and control structure detailed by Table II. For the purposes of the AOAC 
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student officers who executed the Janus simulation, this fictitious enemy force was 
designated as the 168th MRR(-). 

Unit Description Major Systems 



2 Direct Fire Maneuver Battalions BMP 
1 Direct Fire Tank Battalion T80 

4 Indirect Fire Artillery Battalions 2S1, 2S5, BM-21 
1 Air Defense Artillery Company SA7, ZSU-23 
1 Engineer Support Company MTK-LI 

1 Logistics Support Company IF TRU 

1 Helicopter Company HIND (6) 



Table II. REDFOR Task Organization for AO AC 97-3 Janus Simulation 
[Ref. 6] 



B. TERRAIN REPRESENTATION 

The Janus simulation uses digitized terrain developed by the United States 
National Imaging and Mapping Agency (NIMA). The terrain selected for this sim- 
ulation is located within the United States Army’s National Training Center (NTC) 
near Barstow, California, covering a 100 km 2 area with actual military grid coordi- 
nate locations defined by the corner points NK800400, NK790400, NK790390, and 
NK800390. The terrain is displayed in Janus similar to that used by military person- 
nel during actual training and combat operations. It includes contour lines (elevation 
information), roads, rivers, vegetation, and urban areas. The visibility issues as- 
sociated with terrain elevation and visual obstructions are realistically portrayed in 
Janus. These modeled attributes will degrade (or enhance) a combatant’s line-of-sight 
(LOS) and movement. The modeled vehicles and equipment must maneuver under 
the conditions represented by the digital terrain and its overlying natural and man- 
made features. For example, only those vehicles with a “swim” capability for crossing 
water obstacles are permitted to do so without employing bridging assets; and, when 
executed, the vehicles will proceed at a reduced movement rate. 
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The environmental conditions and time of battle are also controlled in Janus. 
The user can specify environmental conditions including temperature, cloud condi- 
tions, and air characteristics affecting smoke propagation. Weather and its effects 
influence the visibility and movement of the modeled combatants. The start time 
can be set in order to execute the combat simulation at a particular time of day. 
This accounts for simulating day and night operations which also affect visibility and 
movement attributes. The environmental conditions used during this simulation were 
summer, desert with 14 kilometers (km) visibility, and a 0800 start time. Since Janus 
accurately models the BLUEFOR and REDFOR weapons systems as a function of 
each system’s capabilities as defined in the database, users must consider all military 
aspects of employing their forces just as they would in actual combat. Invariably, 
Janus will reveal the weaknesses associated with poor planning and reward those who 
have tactically sound plans. As the combat scenario for this research was executed 
by students of the mounted warfare profession, the author assumes the employment 
of tactically sound techniques, procedures, and plans. 

C. SCENARIO EXECUTION SUMMARY 

The initial dispositions of BLUEFOR and REDFOR are depicted in Ap- 
pendix A, Figures 1 and 2 on pages 65 and 67, respectively. The BLUEFOR mission 
is to attack in zone to destroy or cause the withdrawal of the 168th MRR(-) secu- 
rity zone forces who are controlling the key terrain in and around grid coordinate 
NK620130 (Avawatz Mountain Passes) [Ref. 6]. The BLUEFOR executes a sup- 
porting attack in the north followed by a main attack in the south. The REDFOR 
conducts a defensive operation with an objective of stopping the penetration of the 
BLUEFOR and defending the significant terrain as mentioned above [Ref. 6]. The 
REDFOR executes a layered defense with two distinct defensive belts supported by 
deliberate minefields, obstacles, and barriers. 
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At the onset, the BLUEFOR is positioned to conduct offensive operations 
using two distinct avenues of approach (AA). For the purpose of this research, the 
northern and southern AAs will be named AA North and AA South, respectively. The 
main effort of the BLUEFOR attacks along AA South with the supporting attack 
along AA North. As a matter of mounted warfare tactics, the BLUEFOR leads 
with reconnaissance assets in order to locate and direct fires upon the REDFOR’s 
defensive assets. Following closely behind are the BLUEFOR’s initial maneuver and 
engineering assets. This battlefield organization will support direct fire engagements 
with the REDFOR upon their detection by the BLUEFOR reconnaissance elements. 
The BLUEFOR engineering elements are positioned forward to facilitate the clearing 
of obstacles and breaching of REDFOR minefields. Next in the BLUEFOR order of 
battle is the bulk of the fighting forces primarily consisting of tanks and armored 
personnel carriers. The indirect fire artillery assets follow in the BLUEFOR combat 
formation. They are positioned in an attempt to be within range of the REDFOR 
while maintaining reasonable security from REDFOR counter-artillery operations. 
Finally, the support oriented vehicles and equipment make up the rear of the attacking 
BLUEFOR. Positioned well to the rear and in support of the BLUEFOR, six A-10 
close air support (CAS) aircraft are available to be employed as a coordinated part 
of the attack. 

The REDFOR, while maintaining a layered defense supported by maneuver 
obstacles and minefields, also maintains a counterattack capability with reserve forces 
located to the rear of their main defensive area. These combatants can be employed 
to reinforce weakened defensive areas or to counter a penetrating attack made by the 
BLUEFOR. As detailed by Table II, the REDFOR is primarily a mechanized force 
consisting of armored personnel carriers. As such, the REDFOR’s tanks are employed 
judiciously and in positions where they can have the most significant effect on the at- 
tacking BLUEFOR. One of the most significant capabilities of the REDFOR is found 
in the number of supporting indirect fire artillery assets. The REDFOR maintains 



18 



both a close and long range indirect fire capability while keeping the majority of its 
artillery assets outside of the maximum range of the BLUEFOR’s artillery units. The 
REDFOR emphasizes the use of its indirect fire systems and their associated radar 
systems to locate and conduct counter-battery fire against the BLUEFOR. As the 
simulation progresses, the REDFOR attempts to halt the advance of the attacking 
BLUEFOR. 

1. Execution Summary: 0 to 30 Minutes 

The BLUEFOR begins offensive operations by employing its leading recon- 
naissance assets in AA North. Shortly thereafter, the BLUEFOR and REDFOR 
make contact with each other and exchange direct fires. At approximately 11 min- 
utes, the REDFOR long range artillery assets fire special munitions placing two area 
minefields each containing 480 anti-armor mines to thwart the advance of the attack- 
ing BLUEFOR in the north. As of yet, the BLUEFOR attacking along AA South 
have not made contact with the REDFOR defensive positions in that area. During 
and immediately following the BLUEFOR’s minefield breaching operations in the 
north, a heated direct fire battle takes place. Both sides employ their indirect fire 
artillery assets as their reconnaissance assets fix each other’s locations on the battle- 
field. With 28 BLUEFOR and 9 REDFOR combatants killed during this 15 minute 
localized battle, the BLUEFOR attacking along AA North are able to overwhelm the 
REDFOR defensive forces and break through the first defensive line in the north. 
Shortly thereafter, at approximately 27 minutes, the REDFOR long range artillery 
assets reposition to alternate firing positions for security purposes. 

See Appendix A, Figure 3 on page 69 for a pictorial representation of the 
battlefield at the beginning of the scenario. Table III below details the force attrition 
at the completion of the 0 to 30 minute time period. See Appendix A, Figure 4 on 
page 71 for a pictorial representation of the battlefield at 30 minutes. 
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Direct Fire 
Indirect Fire 
Minefields 



BLUE FOR Losses RED FOR Losses 



34 10 

7 4 

4 0 



Cumulative Losses at 30 Minutes 45 14 



Table III. Force Losses for AOAC 97-3 Janus Simulation at 30 Minutes 

2. Execution Summary: 31 to 60 Minutes 

The BLUEFOR, who met with success in the north, proceed east-northeast 
on a high speed AA and are met with little resistance from the REDFOR. At ap- 
proximately 40 minutes, the attacking BLUEFOR moving along AA South come in 
contact with the first layer of REDFOR defenders in the south. The majority of 
the BLUEFOR following the two attacks begin movement north-northeast to take 
advantage of the successful penetration of the REDFOR’s northern defenses. The 
BLUEFOR in contact in the south encounter and begin breaching operations on 
three linear minefields positioned in front of the REDFOR. As BLUEFOR recon- 
naissance forces identify REDFOR targets to their front they direct indirect artillery 
fires against them resulting in only marginal success. At approximately 45 minutes, 
the remainder of the REDFOR’s defending forces in the first defensive layer engage 
the attacking BLUEFOR from AA South. The REDFOR defense in the south is 
much more formidable than that which was found in the north. The BLUEFOR has 
some difficulty identifying penetration points through the REDFOR’s defenses. The 
BLUEFOR suffers a large number of casualties during this battle. The main body of 
the BLUEFOR compresses into a more compact movement formation and increases 
its movement rate. 

Table IV details the force attrition at the completion of the 31 to 60 minute 
time period. See Appendix A, Figure 5 on page 73 for a pictorial representation of 
the battlefield at 60 minutes. 
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Direct Fire 
Indirect Fire 
Minefields 



BLUEFOR Losses REDFOR Losses 



51 14 

9 4 

4 0 



Cumulative Losses at 60 Minutes 64 18 



Table IV. Force Losses for AOAC 97-3 Janus Simulation at 60 Minutes 

3. Execution Summary: 61 to 90 Minutes 

By 65 minutes, the BLUEFOR attacking in AA South are massed at the for- 
ward edge of the linear minefields and conducting passage through the cleared lanes. 
At the same time, six A-10 CAS aircraft enter an aerial holding area for possible 
employment against the REDFOR in the first defensive belt. The REDFOR has a 
small ambush or counter-attack force made up of two tanks and two armored person- 
nel carriers positioned at the extreme south-southeast of its first defensive layer in 
order to block the advancing BLUEFOR from using a high speed A A. However, the 
overwhelming BLUEFOR combat power moving through this area quickly engages 
and destroys this team. At this point, the BLUEFOR begins to reorient its forward 
movement to the north to take advantage of the previous success there. At approxi- 
mately 83 minutes, the BLUEFOR attacking to the north-northeast along AA North 
make contact with the second layer of REDFOR defenses. The BLUEFOR begins an 
indirect artillery barrage against these REDFOR. 

Table V details the force attrition at the completion of the 61 to 90 minute 
time period. See Appendix A, Figure 6 on page 75 for a pictorial representation of 
the battlefield at 90 minutes. 

BLUEFOR Losses REDFOR Losses 
Direct Fire 65 25 

Indirect Fire 10 5 

Minefields 7 0 

Cumulative Losses at 90 Minutes 82 30 

Table V. Force Losses for AOAC 97-3 Janus Simulation at 90 Minutes 
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4. Execution Summary: 91 to 120 Minutes 

At approximately 93 minutes, the six A-10 CAS aircraft depart the aerial hold- 
ing area to attack the second layer of the REDFOR defenses from the southeast to 
the northwest. Shortly thereafter, at 100 minutes, two REDFOR radar controlled 
anti-aircraft systems engage and destroy three A-10 CAS aircraft. The remaining 
aircraft adjust their course and continue their attack on the defending REDFOR 
combatants. Prior to their egress, the remaining three A-10 CAS aircraft are de- 
stroyed by REDFOR shoulder-fired anti-aircraft missiles located in the northeast. By 
105 minutes, the main BLUEFOR attacking elements clear the linear minefields and 
being their movement north-northeast along the high speed AA. At 116 minutes, these 
attacking BLUEFOR make contact and exchange initial fires with the remainder of 
the REDFOR’s second layer defensive forces in the central and southern portion of 
the maneuver area. Heavy exchanges of artillery occur during the development of 
this high intensity battle. 

Table VI details the force attrition at the completion of the 91 to 120 minute 
time period. See Appendix A, Figure 7 on page 77 for a pictorial representation of 
the battlefield at 120 minutes. 

BLUEFOR Losses REDFOR Losses 



Direct Fire 


92 


35 


Indirect Fire 


24 


9 


Minefields 


12 


0 


Cumulative Losses at 120 Minutes 


128 


44 



Table VI. Force Losses for AOAC 97-3 Janus Simulation at 120 Minutes 



5. Execution Summary: 121 to 150 Minutes 

During the battle within the second layer of the REDFOR defense, six ad- 
ditional A-10 CAS aircraft are called forward from their staging area to attack the 
defending forces. The majority of the six aircraft are highly successful until their 
destruction at 131 minutes by radar controlled anti-aircraft fires. Heavy exchanges of 
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direct fires continue to occur and the attacking BLUEFOR divides its formation into 
two distinct attacks, one to the north-northwest and another along the high speed 
AA to the east-northeast. 

Table VII details the force attrition at the completion of the 121 to 150 minute 
time period. See Appendix A, Figure 8 on page 79 for a pictorial representation of 
the battlefield at 150 minutes. 

BLUEFOR Losses REDFOR Losses 
Direct Fire 139 49 

Indirect Fire 32 9 

Minefields 22 0 

Cumulative Losses at 150 Minutes 193 58 

Table VII. Force Losses for AOAC 97-3 Janus Simulation at 150 Minutes 

6. Execution Summary: 151 to 180 Minutes 

The BLUEFOR attack to the north-northwest meets with a tank heavy REDFOR 
defensive element including some reserve forces coming from behind the second de- 
fensive layer. The BLUEFOR reacts by passing additional tank assets through to 
the front of its formation in an attempt to counter the REDFOR tank threat with 
like firepower. At 158 minutes, the BLUEFOR in the northwest directs indirect ar- 
tillery fires against the REDFOR with little effect. At 165 minutes, the main body 
of the BLUEFOR which was following the lead attacking forces moves farther for- 
ward to concentrate fires. The direct fire battle in the north-northwest becomes 
a battle of attrition with the BLUEFOR maintaining numerical superiority. How- 
ever, this translated into a target rich environment for the REDFOR. A successful 
penetration of the defense occurs and BLUEFOR continues their northwestern move- 
ment past the REDFOR’s second layer of defense. As the breakthrough occurs, the 
REDFOR in the south take advantage of the flanking shot opportunities against the 
passing BLUEFOR and are rewarded with a large number of kills. This time period 
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is marked by the BLUEFOR’s largest number of direct fire losses as compared to any 
other 30 minute time period. 

Table VIII details the force attrition at the completion of the 151 to 180 minute 
time period. See Appendix A, Figure 9 on page 81 for a pictorial representation of 
the battlefield at 180 minutes. 

BLUEFOR Losses REDFOR Losses 



Direct Fire 


254 


63 


Indirect Fire 


44 


9 


Minefields 


27 


0 


Cumulative Losses at 180 Minutes 


325 


72 



Table VIII. Force Losses for AO AC 97-3 Janus Simulation at 180 Minutes 

7. Execution Summary: 181 to 211 Minutes 

The BLUEFOR divides into two relatively equally sized formations with one 
attacking north-northeast and the second, south of the first, attacking along the high 
speed AA to the northeast. At 188 minutes, the REDFOR commits its reserves in 
the north to counter the BLUEFOR attack in that sector. After some success, the 
numerically superior BLUEFOR destroys the employed reserves. At 198 minutes, the 
BLUEFOR in the south completely clears the second defensive layer and all of its 
associated minefields. At 206 minutes, the REDFOR commits addition reserve forces 
against the BLUEFOR attacking in the south along the high speed AA in an attempt 
to thwart the advance. At 211 minutes, the combat simulation is terminated simulat- 
ing the achievement of one or more of the BLUEFOR commander’s ending conditions. 
Battle termination standards may include attrition to a specified force level, attrition 
to a specified number of remaining maneuver systems, and/or the accomplishment of 
stated objectives. This particular combat simulation was undoubtedly stopped due to 
a self-imposed time constraint as the simulation was executed by and for the benefit 
of AOAC Class 97-3. 
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Table IX details the force attrition at the completion of the simulation. See 
Appendix A, Figure 10 on page 83 for a pictorial representation of the battlefield at 
the completion of the simulation. 



BLUEFOR Losses REDFOR Losses 

Direct Fire 270 79 

Indirect Fire 45 9 

Minefields 29 0 



Cumulative Losses at Completion Minutes 344 88 

Table IX. Force Losses for AOAC 97-3 Janus Simulation at Completion 



D. UTILIZATION OF THE SIMULATION OUTPUT 

One of the most advantageous properties of the Janus simulation software is 
its ability to capture information about the execution of the combat scenario. The 
force loss totals detailed in the tables above are just samples of the information 
available to the Janus analyst. The Janus post execution output files which are 
important to the attrition rate analysis of this research include the direct fire report, 
coroner’s report, killer/victim scoreboard, and the artillery impacts report. While 
several other reports are available, these four reports will provide the time series of 
casualties data required for the development of attrition rate estimates applicable to 
the ATCAL Phase I algorithm. With the required simulation output data in hand, 
what remains to be developed are the mathematical constructs which serve as the 
basis for the ATCAL Phase I and II algorithms. 
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III. ATTRITION METHODOLOGIES FOR 
AGGREGATED MODELS 

A. BACKGROUND 

Prior to any detailed analysis of the attrition rate estimate process as it applies 
to the ATCAL Phase I and II algorithms (the primary goal of this research), it is 
important to establish a foundation with respect to the various methods for calculat- 
ing attrition and determining the outcomes of battles within aggregated simulations. 
In general, nearly all attrition computation methodologies can be categorized by two 
techniques: the Anderson (Force-Ratio [FR]) model and some type of Lanchester- 
type model. Integral to the second attrition approach are the Bonder-Ferrell method 
and the Com bat Analysis (COMAN) maximum likelihood estimation (MLE) method 
for determining attrition coefficients for use in Lanchester-type equations (of which 
variants are used in the ATCAL algorithms). Although the COMAN MLE ap- 
proach is not specifically used within any current models, its applicability to the 
attrition modeling cannot be overstated. Each of these distinct attrition methodolo- 
gies maintains a relevance to modern aggregated modeling. For example (although 
not all-inclusive), the Anderson (FR) method is one of the two options for assessing 
attrition in TACWAR 4.1, the Bonder-Ferrell method is employed within the United 
States Army’s Vector In Commander (VIC) model, and the ATCAL methodology is 
used in CEM and ITEM. (TACWAR also has an option to use the ATCAL attrition 
methodology rather than FR.) One of the most interesting facts of this research is 
that ATCAL uses none of the attrition methodologies (explicitly) as they are de- 
scribed below. Rather, ATCAL uses an ad hoc method for attrition rate estimation. 
This is truly one of the mysteries of ATCAL for which this research attempts to 
establish the vehicle for further analysis. 
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B. THE ANDERSON (FR) BASED ATTRITION MODEL 

In this theory, a heterogeneous force (its overall strength) can be described 
by one scalar value representing the consolidated individual systems’ contributions 
to the force as a whole. Each system maintains an assigned value (score) relative 
to some common weighting criteria. Summing over all systems multiplied by their 
respective value (score) results in the generation of an overall firepower index (FPI) 
for the force being represented. As this FPI is a single scalar value, the resulting 
process is characterized as a homogeneous force ratio attrition model. The explicit 
representation of the heterogeneous force is replaced by this single scalar measure 
of combat power. This strategy can be expressed mathematically for a hypothetical 
force, X , as shown below in Equation III- 1 . 



FPI* = (Firepower Score) (III.l) 

t=l 

i = integer index each representing a unique weapon system type, total of m 
Sf = value (score) of one weapon system of type i 
X{ = number of weapon system of type i 



Supposing that the FPI of two opposing forces ( X and Y ) are represented by 

FPI* 

FPI^ and FPIy, respectively, then the force ratio, FR^y = , gives a measure of 

FPIy 

relative combat power at the instant at which the values are calculated. The FR^y 
may be used to determine force missions, unit postures, casualty computations, and 
mission success. Specifically, the FR^y can be indexed against personnel loss curves 
to determine attrition to modeled forces. The loss curves employed in this process 
are generally derived from the results of historical battles. This fact is considered to 
be one the inherent strengths of the Anderson (FR) method. On the other hand, one 
of the perceived weaknesses of this method is associated with the determination of 
the individual weapon system values (scores) as discussed below. 
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1. Determination of Weapon System Value 

There is no general agreement in the modeling community on how to deter- 
mine the individual weapon system values (scores), Sf . The fundamental methods 
for determining firepower score values include analysts’ perception of combat value 
(the naive estimate), historical combat performance, measures of weapon firepower 
and support platform characteristics, and, finally, measures of what a weapon can 
kill. This last method attempts to compute the value (score) for a weapon as being 
proportional to the rate at which it destroys the value of opposing enemy systems, 
and, as such, becomes a circularly defined problem. Moreover, this approach’s de- 
pendence on the “rate” at which weapons destroy opposing weapons in battle infers 
the reliance upon assumed or derived non-negative scalar attrition rates (one for each 
weapon system firing at every opposing weapon system for both forces) which account 
for the highly detailed engagement process described in Chapter I. The mathematics 
of this method, called the AntiPotential-Potential (APP) method, amounts to solv- 
ing for the individual weapon system scores through a linear system of equations (in 
matrix form) as detailed below by Equations III. 2 and III. 3. The solutions to these 
equations are found by solving an equivalent system of eigenvalue problems [Ref. 1], 
for which the development is beyond the scope of this discussion. However, as a prac- 
tical example, Appendix B details the application of the Anderson (FR) method and 
the subordinate APP process to hypothetical high resolution output data representa- 
tive of the AOAC 97-3 Janus simulation described in Chapter II. The calculations are 
sampled from the heterogeneous battle between selected BLUEFOR and REDFOR 
weapon systems. Note that this computational excursion does not account for every 
weapon system modeled in the Janus Simulation in an attempt to keep the associated 
eigenvalue problem numerically solvable in closed form without the aid of automated 
software. 
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(III-2) 



CxSf 

CrSj 



E K,,S] 

i=i 

m 

i=l 



(HI-3) 



Force X 
Force Y 

S? 

SJ 

Ktj 

Lji 

Cx 

Cy 



{xi, X2, X3, . . . , x TO } m distinct weapon systems 

{yi 1 V2, yz, ■ ■ ■ , J/n} n distinct weapon systems 

value of one type i system in X force 

value of one type j system in Y force 

matrix of rates at which one Xi system kills ^systems 

matrix of rates at which one Yj system kills ^systems 

proportionality constant 

proportionality constant 



2. Advantages and Disadvantages 

The use of historical data, both in the generation of personnel and equipment 
loss curves and in the determination of weapon system values, generates a certain level 
of model believability and personal comfortableness with the Anderson (FR) method 
of attrition. Moreover, aggregated simulations using this method are credited with 
reproducing comparable results to the recorded outcomes of actual battles. This 
trait is often considered important when attempting to validate and verify a model’s 
usefulness. The desire to extend these result to the modeling of future battles is a 
task which generates a considerable amount of discussion. With respect to determin- 
ing weapon system values (scores), the APP method and its associated eigenvalue 
problem are an attractive procedure for calculating homogeneous force values and 
resulting force ratios for use in aggregated models because of their mathematical ease 
of computation (regardless of their appropriateness). 
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Conversely, there are a significant number of limitations inherent to this pro- 
cedure. A major criticism of this method is that the APP scoring method may be 
inconsistent with the personnel loss curves which are generally based upon historical 
data. Additionally, since Equation III. 1 is additive across weapon systems, there is 
no ability to represent the interactions among complementary weapons on the battle- 
field. This is an especially serious shortcoming considering the nature of the current 
battlefield and the synergies of modern weapon systems. One must exercise caution 
when conducting comparative assessments of force mix or force balance since the 
scalars FPI* and FPIy do not discriminate between balanced and unbalanced forces 
(in terms of size). An FPL* = 100 and an FPIy = 100 are equal forces using this 
method regardless of the number of systems represented or their respective values. 
The reader will also note that Equation III. 1 is linear in the number of weapons Xi 
of each weapon type i. Therefore, a model employing the Anderson (FR) method 
cannot be expected to effectively answer specific questions about the contributions of 
individual weapons. For example, the ramifications of an unmanageably large num- 
ber of a specific weapon may be intuitively obvious, but, the FR model will linearly 
increase the overall FPI without regard to practicality. It is the opinion of the author 
that the overwhelmingly mysterious characteristic of the Anderson (FR) method lies 
within the development of the individual weapon system values (scores). The reader 
may have surmised that in order to apply the Anderson (FR) method, the analyst is 
forced to arbitrarily assume values (scores) for each of the modeled weapon systems 
or to rely upon the killer/victim results of high resolution simulation to develop the 
attrition rates necessary to employ the APP method of solving for weapon system 
values (scores). Associated with using the high resolution simulation results to derive 
the kill rates within the APP method for increased believability, the reader should be 
aware of some of the inherent weaknesses which plague the APP eigenvalue problem. 
The score values derived from the eigenvalue problem are extremely dependent upon 
scenario characteristics and the number of weapons systems which translate into tar- 
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get engagement opportunities. Changes to the attrition rates (the key elements to 
the APP method) will typically cause all weapon system score values to change for 
both forces in ways that are hard to describe or justify. Consider another author’s 
evaluation of the APP method of determining weapon system values: 

The eigenvector score values sometimes change in ways that are hard 
to explain and that have been called paradoxical by some. For example, a 
shift in fire distribution designed to increase the kill rates of a higher value 
enemy target can sometimes reduce the total value index of the firing force. 
Whether this is, in fact, paradoxical depends on how deeply the relationships 
are followed. If there are few of the high value target, then shifting fire away 
from a lower value but more numerous enemy system might very well result 
in a lower total value being killed by the firing force. 

Other anomalies, however, are harder to explain. The numeric values 
of the scores are sometimes oversensitive to small changes in the input kill rate 
matrices... [Ref. 1] 

It is with this knowledge that the analyst is drawn to developing similar 
methodologies using nonlinear value equations, a formidable excursion which is be- 
yond the scope of this research. 

3. Typical Implementation of the Anderson (FR) 
Method 

A generic aggregated combat model, like TACWAR may employ the Ander- 
son (FR) method for attrition calculations in the following way. Assume that a hypo- 
thetical model represents opposing forces down to the division level of organization. 
The model would compute attrition for the divisions in a combat sector with an up- 
date cycle equal to the overall time cycle of the model, possible every 12 or 24 hours. 
Prior to the calculation of attrition, the simulation schedules several events to update 
the forces. These tasks may include movement of forces into and out of the combat 
sectors, allocation of reinforcements, determination of which force is the attacker for 
the cycle, determination of the defender’s posture for the cycle, and computation of 
the results of any scheduled air battles for the cycle. At a point during the execu- 
tion of a generic aggregated combat simulation where the conditions are met for two 
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opposing divisions to engage in battle, the Anderson (FR) method would be invoked 
to determine the respective forces’ FPIs. In a simplified sense, the resultant FRxy 
determines the outcome of the battle with respect to losses. After computing (mostly 
estimating) the parameters required for the Anderson (FR) and subordinate APP 
calculations, the resulting FRxy values are indexed against the appropriate person- 
nel casualty curve to determine overall attrition levels. (Unfortunately, indexing the 
FRxy against the casualty rate and equipment loss curves may not be consistent 
due to differing computational methods. The FRxy is derived using a linear model 
including the APP process, while the casualty rate and equipment loss curves may 
have been derived nonlinearly from historical data.) During each time cycle of the 
simulation (possibly 12 or 24 hours), the APP eigenvalue problem must be recursively 
solved to determine the updated weapon system scores and the resulting forces’ FPIs 
due to their dependence on the number of individual weapon systems remaining. This 
entire process is simply repeated as the simulation continues through additional time 
cycles. 

C. LANCHESTER-TYPE DIFFERENTIAL ATTRITION 
MODELS 

Prior to discussing the Bonder-Ferrell and COMAN MLE methods for esti- 
mating attrition rates, the reader must gain a basic understanding of the overlying 
framework to which the attrition rates are applied. Lanchester-type models origi- 
nated with the early works of F.W. Lanchester, a British engineer and inventor, in 
an attempt to quantify the advantages of concentrating fires in modern warfare. He 
hypothesized that ancient warfare was merely a collection of one-on-one duels where 
force size primarily determined battle outcomes. The resulting mathematical model 
representing this ancient form of warfare is referred to as Lanchester’s Linear Law 
or Firer-Target | Firer-Target (FT|FT) attrition. On the other hand, Lanchester 
supposed that modern warfare maintained the ability to concentrate the effects of 
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multiple weapons on surviving targets, thus creating a many-on-one condition. In a 
hypothetical battle between forces X and Y, this condition would make the force X 
casualty rate proportional to the number of surviving firing systems in force Y. And, 
the resulting casualty-exchange ratio would be inversely proportional to the force 
ratio. This mathematical model for modern warfare conditions is termed Lanch- 
ester’s Square Law or Firer | Firer (F|F) attrition. These underlying Lanchester 
models, with various extensions and improvements, supply the basic framework to 
the ATCAL Phase I and II algorithms, which will be discussed in Chapter IV. 



1. Lanchester’s Linear Law (FT|FT Attrition) 

The primary assumption of ancient warfare is that any weapon system (soldier) 
not engaged in combat must wait until an enemy soldier becomes available before 
joining the battle. The assumption implies that the number of soldiers of force X put 
forth on the battlefield to engage force Y is determined by the number of soldiers of 
force Y present. Therefore, under the conditions of ancient warfare, there should be no 
advantage to concentrating forces. A particular soldier from a force is either engaged 
in a one-on-one duel, or not. Intuitively, each force size is a function of time, where 
the continuous real time variables x(t), y(t ), and t are approximations to their discrete 
analogs in an actual battle. The change of each force size with respect to time depends 
upon the number of firers engaged in combat and can be represented mathematically 
by the (homogeneous) Lanchester Linear Law (FT|FT), Equations III. 4 and III.5, 
where and represent the change in the size of forces X and Y, respectively, 
over time. 



dx 

— = - ax(t)y(t ) 


with ar(0) = Xo 


(III.4) 


^ = -bx(t)y(t) 


with y(0) = 2 /o 


(IH.5) 


x(t ) = size of force X at time t 
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y(t) = size of force Y at time t 

a = attrition rate at which an individual Y firer 
kills X targets 

b = attrition rate at which an individual X firer 
kills Y targets 



2. Lanchester’s Square Law (F|F Attrition) 

Extending the assumptions of ancient warfare to modern battles, one gains the 

attribute that many firers can engage a single target, thus achieving the advantages 

associated with concentration of fires. Two general cases of combat developed out 

of this study, namely, area fire and aimed fire. First, Lan Chester assumed area fire 

to be a situation where a force spreads its fires over a general area occupied by the 

enemy rather that aiming directly at individual enemy weapon systems. Assuming 

that these area fires are uniformly distributed over the target area and that the 

area is independent of the number of targets, the Lanchester Linear Law (FT|FT) 

was assumed to hold true for area fires. Second, aimed fire assumes that individual 

targets are identified and attacked by one or more opposing firers. In this case, 

Lanchester assumed that the attrition of force X depended on how many force Y 

firers were directing aimed fires at force X. Mathematically, the basic homogeneous 

form of Lanchester’s Square Law (F|F) can be described by a system of differential 

dx dy 

equations, Equations III. 6 and III. 7, where — and — represent the change in the 

at at 

size of forces X and Y, respectively, over time. 



dx 

dt 


-ay(t) 


with x(0) = x 0 


(111-6) 


II 


—bx(t) 


with y(0) = y 0 


(III.7) 


x{t) = 


size of force X at time t 




y(*) = 


size of force Y at time t 
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a = attrition rate at which an individual Y firer 
kills X targets 

b = attrition rate at which an individual X firer 
kills Y targets 

The homogeneous representation of Lanchester’s Square Law (F|F) can be ex- 
tended to meet the needs of heterogeneous forces by accounting for all of the possible 
firer/target combinations in the simulated battle. The resulting heterogeneous for- 
mulation of Lanchester’s Square Law (F|F) is provided by the system of differential 
equations, Equations III. 8 and III. 9, shown below. 



dr- n 

-TT = “ ( IIL8 ) 

^ j= 1 

dv 171 

-%■ = -Em '**(<) ( In - 9 ) 

ai 2 — 1 



Vid) 

(Lij 

bji 
vl( .. 

X force 
Y force 



size of force X, system i, at time t 

size of force Y, system j, at time t 

attrition rate at which Yj firers kill X{ targets 

fraction of Yj firers allocated against X t targets 

attrition rate at which X{ firers kill Yj targets 

fraction of X{ firers allocated against Yj targets 

{x u x 2 ,x 3 ,...,x m } 

{yi,y2,y3,---,y n } 



The reader should be aware that Equations III. 8 and III. 9 represent one pos- 
sible form of the Lan Chester Square Law (F|F) equations. While the equations pre- 
sented represent work which was primarily accomplished during the 1960s, the overall 
methodology has evolved over time. In the 1970s, fire allocation factors which influ- 
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ence the attrition rates were developed by submodels executed in support of the 
overall process. Nevertheless, the equations described above appear in a fundamental 
form which aids in the understanding of this topic. 

3. The Significance of a,ij and bji 

As Lanchester’s original purpose was a qualitative justification of the principle 
of concentration as an advantageous tactic in modern warfare, the derivation of the 
heterogeneous attrition rate coefficients a,j and bji was of no particular concern. How- 
ever, today, as Lanchester based attrition models are employed as part of theater level 
aggregated models, the significance of these attrition rates cannot be overemphasized. 

As part of the model verification process, current aggregated models must provide 
reasonable and quantitative results. These numerical values assigned to a^- and bji 
are what represent the synergies of battle between various weapons systems and the 
very nature of specific actions on the battlefield. Their derivation is what brings 
believability to the equations present within combat models. Several methods have 
been proposed and utilized over the years to derive the attrition coefficients. The 
remaining background discussion of this chapter focuses on two of these methods of 
attrition rate estimation which are currently in use. 

D. THE BONDER-FERRELL METHOD FOR ATTRITION 
RATE ESTIMATION 

The Bonder-Ferrell method for estimating attrition rates makes use of an in- 
dependent analytical model which does not depend on input from other sources. 
Instead, a stochastic process model is developed for each Yj firing at a single passive 
type Xi target, and the expected time between casualties (E[T-j]) is determined, where 
T-j is a random variable representing the time it takes a single Yj firing system to 
kill one Xi target. Assuming that the casualty process is a Markov (Poisson) process 

(see [Ref. 7]), the attrition rate parameters can then be expressed as a,y = 3, y and 

E i T ij\ 
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bji = . While Bonder’s original work focused on deriving the equation to com- 

^ l* jii 

pute E[T] by calculating the distribution of the number of shots required to achieve 
a specified number of hits, Ferrell’s extension of this work modeled this process as a 
Markov renewal process (see [Ref. 7]) with three discrete states (new engagement, hit, 
and miss) and the associated transition probabilities and times. Thus, the expected 
time required to kill an enemy target is the long run expected time required within 
the Markov renewal process until moving on to a new engagement. As early studies 
in this area focussed on homogeneous forces, continued work in the 1970s extended 
the concept to heterogeneous battles by combining conditional kill rates with target 
priority lists. Additionally, target acquisition was characterised as either parallel or 
serial. 



One of the most significant advantages to the Bonder-Ferrell method of at- 
trition coefficient estimation is that the computations are based upon what is going 
on in the battle at the given moment. There is no requirement to reference a data 
base library of high resolution simulation output to match up the conditions of battle 
between the overlying model and the high resolution battle. Additionally, since the 
underlying data used in these computations are basically measurable, the results are 
more credible to the user. This method is generally considered to be restrictive in 
that the in-depth engagement model is analytic, thus suppressing some detail. This 
form of attrition rate estimation also precludes the ability to capture synergistic ef- 
fects of weapon system combinations as the underlying process does not account for 
weapon system interactions. The reliance upon engineering level data does lend this 
approach to the explicit definitions of engagements; however, the amount of data 
required on each weapon system is great. Relying upon closed form equations for the 
development of the a,j and bji , the Bonder-Ferrell method allows for the re-evaluation 
of the attrition coefficients at each time step of the simulation by simply recalling and 
recomputing the required input values. 
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E. THE COMAN MAXIMUM LIKELIHOOD ESTIMA- 
TION (MLE) APPROACH TO ATTRITION RATE 
ESTIMATION 

The COMAN MLE attrition estimation methodology is a fitted parameter 
model which examines the time series of casualty times (from a high resolution sim- 
ulation) and computes the maximum likelihood estimate for the mean time between 
casualties. While other statistical methods exist for determining point estimates of 
this type (including the method of moments, Bayesian estimation, and the method of 
least squares), the MLE approach is the only method which has had any significant 
application in combat simulations [Ref. 8]. Similar to the Bonder-Ferrell approach, 
the COMAN MLE method relies upon the assumption that the time between casu- 
alties is a Markov (Poisson) process with the memoryless property (time until the 
next casualty does not depend upon the previous time before the last casualty). The 
COMAN MLE method is most easily explained by detailing its mathematical con- 
struct. The underlying assumption of this method is that the high resolution simula- 
tion data represents a sample from the continuous analogue to the deterministic homo- 
geneous equations of Lanchester’s Square Law (F|F) (see Equations III.6 and III. 7), 
i.e., the continuous-parameter Markov chain with transition probabilities given by 
Equations III.10 and III. 11 [Ref. 9]. 



Prob [A casualty in a small time interval of length At] = an(t)At 

(III.10) 

Prob[E casualty in a small time interval of length At] = 6m(t)At 

(mu) 



n(f) = the number of Y combatants at time t 
m(f) = the number of X combatants at time f 
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Additionally, consider the following random variable definitions and their as- 
sociated properties required for the development of the MLE attrition coefficient es- 
timation process. The reader can assume that any realization of a described random 
variables will be notationally represented by its corresponding lower-case letter. 



c? 



c? 



r x 

°T 



4 



4 + 4 



1 if the k th casualty is an X combatant, 

< 

0 otherwise, 

1 if the k th casualty is an Y combatant, 

< 

0 otherwise, 

K 

y t c x total number of X casualties 

k= i 

K 

y, 4 total number of Y casualties 

k—l 

K total number of X and Y casualties 



n(tk) 

m(t k ) 



n k the realization of the number of Y combatants 
just after the occurrence of the k th casualty 

rn k the realization of the number of X combatants 
just after the occurrence of the k th casualty 



Interpreting the random variable definitions given above, note that during 
the time interval there are m k X combatants and n k Y combatants alive 

for values of k = (0, 1, 2, . . . , K — 1). Armed with high resolution simulation output 
data of the form {(tj, c x , 4 ), fa, 4 ■> 4) ■> • • • , i^K, c /f)}, the statistical estimates 
for the continuous time Markov chain representation of the deterministic form of 
Lanchester’s Square Law (F|F) can be computed using the method of maximum 
likelihood estimation (MLE). As a battle develops, the casualty processes for the X 
and Y forces are simply two superimposed Poisson processes with respect to time (a 
random variable). Developing and maximizing these likelihood functions from their 
applicable joint probability density function (PDF) (representing the times between 
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casualties) will result in the derivation of the attrition rate estimates, a and b. The 
maximum likelihood estimation process can be summarized by the steps detailed in 
Table X. 



Step 1 


Find the probability density function for the time until an X casualty. 
Find the probability density function for the time until a Y casualty. 


Step 2 


Construct the respective likelihood functions. 

(i.e., the density functions for the observed sequences of events) 


Step 3 


Determine the values of a and b that maximize the 
respective likelihood functions. 



Table X. Maximum Likelihood Estimation Process Steps 

[Ref. 9] 



1. Maximum Likelihood Estimation Steps 

Applying the three maximum likelihood estimation steps to the stochastic rep- 
resentation of Lan Chester’s Square Law (F|F) represented by Equations III. 10 and III. 1 1 
will result in the statistical derivation of a and b. The author has chosen to detail 
the development of the estimated attrition parameters a and b from the X casualty 
point of view. Initially, C* and 1 — C * will be used to represent X and Y casualty 
occurrences, respectively. 

a. Step 1 

Let the random variable Sx denote the time between X casualties with 
the well-known exponential PDF fs x (s x ) = an exp _ ^ an ^ x . Similarly, the PDF for 
the time between Y casualties is fs Y (s y ) = bmex p - ( 6m W. Further, let the random 
variable S denote the time between any two consecutive casualties be defined as the 
minimum of Sx and Sy. The resulting PDF describing the time between any two 
casualties is 

fs(s) = (an + bm) exp- (an+6m)s . (III.12) 



41 



Incorporating the previously defined variable C * , which equals one when the k th 
casualty is an X combatant and zero when the k th casualty is a Y combatant, one 



can further specify Equation III. 12 above. Note that C* = 1 if Sx < Sy. The 



probability of this occurrence is Prob^jf = 1) = 



an 



an + bm 



Likewise, a Y casualty 



occurs if Sx > Sy which implies Cjf = 0 and Prob(C/ = 0) = 
the joint probability of the random variables S and C* is 



bm 



an + bm 



Therefore, 



Prob(5' <s,C? = i ) = (1 - exp-( an+6m ) s ) • Prob(c£ = i) 



where 0 < s < oo and i = 0, 1 



= ( 



an 



an + bm 



r •( 



bm 

an + bm 



i-ci 



(1 - exp- (an+!>m)s ). 



Therefore, the representative joint PDF of S and C jf is 



fsc x ( s i c k ) = ( — ~T — Y* ' ( — ~~V — )* c * • (an + bm) 
j$,c k \,kj K an + bm J K an + bm J K ’ 

= (an) c * {bm) l ~ c * exp _ ^ an+fcm ^ 



exp- (an+6m)s 



b. Step 2 

The likelihood function is developed by considering that the K casu- 
alties have occurred at times (ti,t 2 ,t 3 , . . . ,tx) with a total number of X and Y ca- 
sualties represented by and c^, respectively, where + Cj = K. Recalling that 
C$ = 1 C% = 0 and C* = 0 = 1, Equation III. 13 can be functionally 

written as 



(arc*_i) c * (6mjb_i) c * exp- (on *- 1+fem ' ! - ,)(t *- t *- l) . (III.14) 
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Equation III. 14 represents the contribution of the k th casualty to the likelihood func- 
tion. As the times between casualties are independent random variables, the likeli- 
hood function for the observed sequence of events, Equation III. 15, is the product of 
all of the independent contributions. 

K 

L(a,b ) = ex P“ (oB ‘" 1+6m ‘" l)(tfc "*‘" ,) 

A:=l 

(III. 15) 



c. Step 3 

Finally, the COMAN MLE process requires the maximization of this 
likelihood function. However, it is a common mathematical practice to maximize the 
natural logarithm (In) of the likelihood function as detailed below by Equation III. 16. 
Both maximum values will occur at the same point [Ref. 10] and the logarithmic form 
is more easily manipulated. 



K K K 

ln[L(a, 6)] = £ c x ln(an*_i) + Y, c l ln(6m*_i) - X^( an fc-i + - t k -i) 

k= 1 k—\ k=l 

(III.16) 

Continuing with Step 3, one can take the partial derivatives of Equation III. 16 with 
respect to a and b to get Equations III.17 and III. 18, respectively. 



flln[£(q,6)] 

da 

d\n[L{a, 6)] 

db 



c x K 
a k = i 

c? K 

-r 

0 k=l 



(111.17) 

(111.18) 
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Finally, maximize Equations III. 17 and III. 18 by setting each of the partial deriva- 
tives equal to zero and solving for a and b. The resulting solutions for a and b, 
denoted a and S, maximize the likelihood function (Equation III. 15) and serve as the 
COMAN MLE estimates for the attrition rate coefficients. See Equations III. 19 and III. 20. 



£r 

£fc=l Tlk-l(tk — h-i) 



J2k=i mk-i(tk tk— i) 



(111.19) 

(111.20) 



2. Summary of the COMAN MLE Method 

While the presented mathematical description of the COMAN MLE method 
for estimating attrition rate coefficients is applicable only to the homogeneous rep- 
resentation of Lanchester’s Square Law (F|F) in one of its simpler forms, the reader 
can see that the general method is quite simplistic. One of the simplifying charac- 
teristics of the developed example is the appearance of the attrition rate coefficients 
as linear terms only. It is quite possible that the true nature of the attrition rate 
coefficients is non-linear in actual combat; however, for the purposes of this dis- 
cussion, the author has chosen to represent the mathematics of the method with 
tractable forms. As a practical example, Appendix C details the application of the 
COMAN MLE method for attrition rate estimation to the high resolution output data 
from the AOAC 97-3 Janus Simulation described in Chapter II. The calculations are 
presented for the mechanized battle between the BLUEFOR M2 (with TOW missile) 
and REDFOR BMP-2 (with AT-5 missile). Note that this computational excursion 
makes use of Lanchester’s Square Law (F|F) (as the mechanized battle is one char- 
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acterized by exchanges of aimed fire) and the mathematical framework leading up to 
and including Equations III. 19 and III. 20. 

F. EXTENSION OF ATTRITION RATE CONCEPTS 
TO THE ATCAL ALGORITHMS 

By no means should the reader consider this treatment of attrition rate anal- 
ysis to be exhaustive. A significant number of noteworthy modifications have been 
proposed and are used within Lanchester-type deterministic models. For a highly 
detailed inspection of these modifications, the interested reader should consult the 
prolific works of Professor James G. Taylor, Naval Postgraduate School, and co- 
advisor to this research. However, with the general background development of the 
methods for attrition rate estimation complete, the reader is now prepared to consider 
their application to (or lack thereof) the ATCAL algorithms. The specific modifica- 
tions to the underlying Lanchester-type model which form the basis for the ATCAL 
methodology can be primarily attributed to Gordon M. Clark and his PhD Thesis, 
“The Combat Analysis Model” [Ref. 8], while ongoing analysis is being conducted by 
James G. Taylor [Ref. 2]. However, recall that the specific attrition rate parameter 
estimation process used in ATCAL Phase I can be described as an ad hoc process at 
best. 
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IV. THE UNDERLYING 
LANCHESTER-TYPE MODEL FOR ATCAL 

As the primary purpose of this research is to establish a test mechanism for 
evaluating the current ATCAL attrition algorithm as it applies to aggregated combat 
models (like CEM and TACWAR), the interested reader must develop a familiarity 
with ATCAL’s underlying Lanchester-type models and their application via explicit 
assessment equations. As detailed documentation describing the theoretical bases of 
ATCAL is not available from its original author, this chapter attempts to outline the 
construct of the ATCAL attrition algorithm based primarily upon a very thorough 
treatment of the subject by James G. Taylor. [Ref. 2] It is with an understanding of 
the current ATCAL algorithm that the interested reader could conduct research for 
improving its mathematical efficiency or developing mathematical variants to fulfill 
other combat model requirements. The reader may find it helpful to review Chap- 
ter III, Section C, which details the basic framework of Lanchester’s Linear (FT|FT) 
and Square (F|F) Laws because of the parallelisms which extend to the Lanchester- 
type models present in the ATCAL attrition algorithm. 

A. THE HOMOGENEOUS FORCE MODEL 

Considering the case of homogeneous combat between two opposing forces, the 
ATCAL attrition algorithm employs the nonlinear Equations IV.l and IV.2 which can 
be attributed to Gordon M. Clark [Ref. 8]. 



dx 

dt 


-0(1 -?*<•>)»(() 


with x(0) = x 0 


(IV.l) 


II 

~<3 |^3 


1 

1 

<ce 

V ?T 

es» 


with 3 /( 0 ) = j/o 


(IV.2) 


H 

V 

o 


size of force X at time t 
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y(t) > 0 

a > 0 

P > 0 

0 < P x{t) < 1 

0 < g y{t) < 1 



size of force Y at time t 

rate at which an individual Y firer kills acquired X targets 
rate at which an individual X firer kills acquired Y targets 
probability that an individual Y firer has no 
X target available for engagement 
probability that an individual X firer has no 
Y target available for engagement 



The reader may find it more convenient to express Equations IV. 1 and IV. 2 
with respect to target availabilities rather than non availabilities. To do so, let 1 — (1 — 
A)*^) = 1 — p x v\ which is the probability that an individual Y firer has one or more 
X target(s) available for engagement. And, analogously, let 1 — (1 — B) y ^ = 1 — q y ^K 
With these substitutions, Equations IV. 1 and IV. 2 can be expressed as follows. 



— = — a[l — (1 — with x(0) = xo (IV.3) 

at 

j t = -p[\ - (1 - B) y U]x(t) with y(0) = yo (IV.4) 

Considering the elements of Equation IV. 3, the reader will note that the mul- 
tiplication of the conditional kill rate (a) with the probability of target availability 
(1 — (1 — produces the expected rate at which an individual Y firer kills X 

targets. Similarly, the multiplication of ft and 1 — (1 — B ) y W yields the expected rate 
at which an individual X firer kills Y targets. 

It has been shown that force levels can never become zero or negative as a 
matter of theory (see [Ref. 9]). However, as the attrition calculations are necessarily 
executed on computer platforms, the potential for achieving “computer zero” must 
be recognized. So, in practice, the attrition models above must be defined to be valid 
only when the forces levels, x(t) and y(t) are greater than zero. This refinement of the 
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homogeneous attrition model within the ATCAL algorithm produces Equations IV. 5 
and IV. 6. 



dx 


— a[l-(l -AfW)y{t) 


for x(t),y(t ) > 0, 


(IV.5) 


dt 


0 


otherwise. 


dy 1 

dt | 


\ -/?[!-(! -£)*>]x(i) 


for x(f),y(t) > 0, 


(IV. 6) 


1 0 


otherwise. 





This positive force level requirement can be easily incorporated in computer algo- 
rithms by including a conditional check of the current force level minus any calcu- 
lated change to the force level. If the resulting force level is less than zero, then the 
algorithm should return a value of zero. Otherwise, the current force level is updated 
by subtracting the calculated attrition. 



1. Application of the Averaging Operator 

In order to transform Equations IV. 5 and IV. 6 into explicit assessment equa- 
tions as employed in the ATCAL algorithm, it is important to develop the mathe- 
matical properties of the averaging operator (over time) given by 

l r* 

Ave[x\ = x = - / x{s)ds (IV. 7) 

t Jo 

where x(t) is a function of time. 



Application of this operator to the Lanchester-type equations developed above will 
result in an operator equation for the force levels, x(f) and y(i), respectively. Define 
the number of casualties inflicted on the X force as Ax and the number of casualties 
inflicted on the Y force as Ay. 



Ax = Xo — x(f), where Xo — x(0) 

Ay - y 0 — y{t), where y 0 = y(0) 

dx 

Now, consider applying the averaging operator to — which results in 

at 



,dx. 1 A dx 

Ave — = - I — ■ ds 

1 dt 1 t Jo ds 
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£ So] 

[s 0 - s(<)] 

t 

Ax 

~ r 

Hence, applying the averaging operator to the homogeneous Lanchester-type 
Equations IV. 5 and IV. 6 representative of the ATCAL algorithm, one obtains the 
operator equations, Equations IV. 8 and IV. 9 below. 

= -o[l - (1 - A)^]y(t) 

= Ave[—a[ 1 — (1 — A) I ^^]t/(i)] 

= —aAve^l — — A) x ^}y{t)\ 

= atAve\[ 1 — (1 — A) x ^]y(2)] (IV.8) 

and, similarly, 

Ay = (3tAve[[\ - {\ - B) y{t) }x{t)\. (IV.9) 

The reader will note that the arguments of the averaging function appearing in 
Equations IV. 8 and IV.9 above are highly nonlinear. Unfortunately, these equations 
cannot be further simplified to explicit expressions involving only average force levels 
without making some type of approximation. The key approximation made within the 
ATCAL algorithm is that the average of this nonlinear argument in x(t ) and y(t ) is 
equal to the same functional relationship between x and y, the results of the averaging 
function applied to x(t) and y(t). Although this approximation is not precise, it is 
often applied to simplify averages of nonlinear functional relationships. It is this 
simplification that allows ATCAL to move forward computationally. Equations IV.8 
and IV.9 can now be replaced by Equations IV. 10 and IV. 11. 



Ave[ 



dx 

dt 
dx . 

dt' 

Ax 



t 

Ax 



Ax = at[\ — (1 — A) x ]y 



(IV.10) 
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(IV.ll) 



Ay = pt[ 1 - (1 - Bf)x 

x = time average of the force level x{t) 
y = time average of the force level y(t) 

Ax = number of X casualties during time t 
Ay = number of Y casualties during time t 

2. Applying the Results for Homogeneous Forces 

As written, Equations IV. 10 and IV.ll have four unknown quantities, namely, 
the average force levels (x and y) and the numbers of casualties (Az and Ay). With 
two equations and four unknowns, additional assumptions are required in order to 
solve this system of assessment equations. The second key assumption made in the 
ATCAL algorithm is that average force levels can be represented by the expressions 
v — a r. — [Ref. 11] It so happens that these expressions 



x = 



a and y = 

H 1-ff) 



ln( 1-^) 

^0 / v 3/0 ' 

hold true, respectively, if x(t) and y(t) are exponential functions (increasing or de- 



creasing) over time. As force levels necessarily decrease during battle, the negative 
exponential function is assumed. Combining this information with Equations IV. 10 
and IV.ll, the underlying Lanchester-type attrition model used to determine the 
numbers of casualties during each assessment period is represented in the ATCAL 
algorithm by the following system of four simultaneous nonlinear equations in the 
four unknowns, x,y, Ax, and Ay. 



Ax 


= at[l — (1 


Ay 


= /?*[!-(! 


X 


— (~ A ;) 
'"(i-ff) 


7 / 




y 


is 

i 

jT 

1 



While small perturbations of this notation are explicitly modeled in the ATCAL 
algorithm, this system of equations should be thought of as the basic ATCAL casualty 
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assessment system of equations. This form better suits future research in applying 
numerical solution techniques to solve the system (recognizing the requirement for 
excluding negative force levels as described in Equations IV. 5 and IV. 6). The actual 
method used in ATCAL to solve this system of nonlinear equations is a heuristic 
form of successive substitutions which is beyond the scope of this research. The 
interested reader may find the study of alternative numerical solution techniques 
extremely relevant to verifying current ATCAL methodologies and proposing future 
replacement ones. 

B. THE HETEROGENEOUS FORCE MODEL 

Considering the case of heterogeneous combat between two opposing forces 
made up of m distinct X force weapon systems and n distinct Y force weapon systems, 
the ATCAL attrition algorithm employs the nonlinear Clark Equations IV. 12 and 
IV. 13 detailed below. In these equations, i = 1 denotes the lowest priority target 
type and i = m denotes the highest priority target type within the X force. And, 
j = 1 denotes the lowest priority target type and j — n denotes the highest priority 
target type within the Y force. Again, the reader may find it helpful to compare 
these equations with Equations III. 8 and III. 9 of Chapter III. 



dxi 

dt 

dyj_ 

dt 

X{(t ) > 0 
yj(t) > 0 
*i(0) 





-p« w ) 


n 


rlfvim 


for i — 


1, 


2,- 


■ • , m 


3 = 1 




k=i + 1 


























(IV. 12) 


m 


- ?r’) 


n 

n 


(t) 


for j = 


1, 


2,- 


••,n 


1 = 1 




k=j + 1 


























(IV. 13) 


number of 


weapon 


systems of type i 


in force X 


at 


time t 


number of 


weapon 


systems of type j 


in force Y 


at 


time t 



= x,o for i = 1 , 2, • • • , m 



52 



w(°) 

Q,j > 0 
0H > 0 
0 < < 1 



0 < $ (t) 



< 1 



Vjo for i= 1,2, ■•■,n 

rate at which an individual Yj firer kills acquired Xi targets 
rate at which an individual Xi firer kills acquired Yj targets 
probability that an individual Yj firer has no 
X{ target available for engagement 
probability that an individual Xi firer has no 
Yj target available for engagement 



Similar to the homogeneous case, the reader may find it more convenient to express 
Equations IV. 12 and IV. 13 in terms of target availabilities by letting 1 — (1 — Aij) x ‘^ = 
1 — p*j^ and 1 — (1 — Bj{) y ’W = 1 — q V jl^ ■ With this substitution, Equations IV. 12 
and IV. 13 can be expressed as follows. 



J T . n m 

ai j = 1 k=i + 1 

(IV. 14) 

j m n 

-f = -EftiK'-O n 0 - i = l, 2 ,--.,n 

al i= 1 k=j+\ 

(IV. 15) 

Again, in order to preclude the possibility of achieving “machine zero” (or 
negative force levels) during numerical attrition calculations, Equations IV. 14 and 
IV. 15 can be further specified so that the force levels X{(t) and yj(t) are valid only 
when they are greater than zero. One final notational simplification renders the 
heterogeneous formulation of the underlying Lanchester-type model for ATCAL more 
compact. Consider the set of indices corresponding to target types with a higher 
priority than a given target type by formally defining 

Iij = {k\k € integers} 

X target type k has a higher priority than target type i for a Yj firer; and, 
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Jji = {A:| A: 6 integers} 

Y target type k has a higher priority than target type j for a Xi firer. 

With this simplifying set notation, a preferred form for the underlying Lanchester- 
type equations for ATCAL appear below in Equations IV. 16 and IV. 17. This rela- 
tively simple representation results in less computational requirements when imple- 
mented as the attrition model of a large-scale aggregated combat model. 



d t x 

dt 



djV 

dt 



' <*<;{(! -pS (,| )n 4s / u p3" , }w(() 

0 



for Xi(t),yj(t) > 0, 
otherwise. 



- zt, A-<{( i - $ t0 ) rw,. 

0 



(IV. 16) 

for Xi(t),yj(t) > 0, 
otherwise. 



(IV. 17) 



1. Application of the Averaging Operator 

Without loss of generality, the application of the averaging operator (with 
associated notational definitions) to the heterogeneous Lanchester-type equations will 
result in operator equations for the force levels, x;(f) and yj(t). In order to maintain 
consistent notation, the author has chosen to represent the equations using the “target 
availability” notation of Equations IV. 14 and IV. 15. 



n (i-^r w tew 



3 - 1 



k=i + 1 



Ax,- 






3 = 1 



k=i + 1 






i=i 



k=i+l 
m 



A Xi = 



(£a 0 /!ve[{(l-(l -*,)"<•>) II (1 ~ 

j= 1 A: = z -h 1 

(IV. 18) 
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and, similarly, 



t=l *=7+1 

(IV. 19) 

Analogous to the homogeneous case, the key approximation made within the 
ATCAL algorithm is that the average of these nonlinear arguments in x,-(£) and yy(£) 
are equal to the same functional relationships between x; and yy , which are the results 
of the averaging function applied to x,(f) and yy(£). Applying this approximation, 
Equations IV. 18 and IV. 19 can be replaced by Equations IV. 20 and IV. 21. 



Ax,- = n (!--+,)**}» (IV .20) 

j—l k=i + 1 

771 71 

A Vi = <Eft,{(l-(l-fy 5 ’) n (IV .21) 

1=1 k=j + 1 

2. Applying the Results for Heterogeneous Forces 

Similar to the homogeneous argument presented in the previous section, both 

the average force levels ( x; and yy)and the numbers of casualties (A Xi and Ayy) 

are unknowns in Equations IV. 20 and IV. 21. Therefore, this system is made up of 

(m + n ) nonlinear equations in 2 • (ra + n) unknowns. In order to increase the number 

of equations to equal the number of unknowns , ATCAL again makes use of an 

(—Ax;) 



assumption for the average force levels resulting in the expressions x,- = 

(-Aw) 



'"(1 - 



and yj = 



ln(l -£&■)' 

v Vjo ' 



[Ref. 11] Combining this information with Equations IV. 20 



and IV.21, the underlying Lanchester-type attrition model used to determine the 
numbers of casualties during each assessment period is represented in the ATCAL 
algorithm as the following system of 2 • (m + n ) simultaneous nonlinear equations in 
the 2 • (m + n) unknowns, x;,yy, Ax;, and Ayy. 
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Ax,- — 1 5^j_i (1 A{j) X ') rijt=t+l(l A kj) Xk }yj 

A Vi = tJ2i = \ f3ji{{] — (1 — rijfc=j+lO — Bki)* k }%i 



While small perturbations of this notation are explicitly modeled in the ATCAL 
algorithm, this system of equations should be thought of as the basic ATCAL caualty 
assessment equations. This form better suits future research in applying numerical 
solution techniques to solve the system (recognizing the requirement for excluding 
negative force levels as described in Equations IV. 16 and IV. 17). The actual method 
used in ATCAL to solve this system of nonlinear equations is a heuristic form of 
successive substitutions which is beyond the scope of this research. The interested 
reader may find the study of alternative numerical solution techniques extremely 
relevant to verifying current ATCAL methodologies and proposing future replacement 
ones. 



C. SUMMARY OF THE LANCHESTER-TYPE ASSESS- 
MENT EQUATIONS IN ATCAL 

The homogeneous and heterogeneous systems of nonlinear equations found 
at the ends of Sections A and B, respectively, comprehensively describe the under- 
lying Lanchester-type models and their basic assessment equations for the ATCAL 
attrition methodology. “The most significant implication of these results, however, 
is that a much more convenient way of generating numerical results for ATCAL 
has immediately been established: namely, numerical integration of the underlying 
Lanchester-type differential equations.” [Ref. 2] Moreover, with these mathematical 
constructs, the interested reader is now prepared to further investigate the soundness 
of the heuristic successive substitution method of ATCAL and compare it to more 
robust numerical solution techniques. 
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V. 



CONCLUSION 



A. CONCLUSIONS 

This research has developed a practical test mechanism for studying the attri- 
tion methodologies of aggregated combat models. The Janus high resolution simula- 
tion software, capable of modeling brigade sized (and below) operations, is completely 
adequate for producing the required high resolution output data required for the anal- 
ysis of aggregated attrition methodologies. Janus is both widely available within the 
military community (to include the Naval Postgraduate School) and is easily learned 
through the use of a tutorial. See [Ref. 5]. This user friendly computer simulation 
provides the researcher with sufficient detail and a mechanism for capturing detailed 
interactions on the modeled battlefield. The ability to capture the highly detailed 
interactions between weapon systems is what qualifies Janus as a valuable tool for 
conducting analysis of the attrition algorithms inherent to ATCAL (and other ag- 
gregated attrition methodologies). The reader should follow the steps outlined in 
Section C to initiate such a study. 

B. RECOMMENDATIONS 

First, having established the vehicle (Janus) for generating the data required 
for developing attrition (kill) rate parameters, the author has opened the door for 
continued research in the area of attrition methodologies. As this research developed 
and discussed one such combat scenario, the continued study of attrition calculations 
would require additional Janus scenarios with varying parameters. The creation and 
execution of such scenarios could easily be incorporated into current high resolution 
modeling courses taught at the Naval Postgraduate School. 

Second, as discussed in Chapter IV, the nonlinear attrition assessment equa- 
tions resulting from the Lanchester-type equations used in the ATCAL algorithm are 
currently “solved” using a heuristic method of successive substitutions. Significant 
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follow-on research is applicable in the employment of a more robust numerical so- 
lution technique for these systems of equations. This area of continued study does 
not necessarily depend upon additional Janus scenarios; however, additional scenarios 
would provide a data source for technique verification. 

Third, the computer code which contains the ATCAL Phase I and II algo- 
rithms exists in several mediums. A significant endeavor integral to making ATCAL 
(and any proposed modifications) more efficient includes a thorough analysis of the 
current computer code which exists in at least two languages, FORTRAN and C. 
The interested reader may find that the algorithms, which make up ATCAL would 
be more efficient in one or more specific computer languages. For example, a proposed 
numerical solution technique may be more conducive to one computer language than 
another. 

Finally, an automated bridge is required between the Janus high resolution 
output data and the computer code used to perform the aggregated attrition calcula- 
tions (like ATCAL). While this capability may not be directly applicable to the actual 
execution of aggregated models as they currently exist, it is integral to the efficient 
study of the algorithms themselves. It is in the research arena that this automated 
capability for linking Janus output to attrition algorithms would be helpful. JETS 
may be an important tool for continued study of this recommendation. 

As combat models are an integral part of the military analysis and training 
communities, continued analysis of attrition methodologies (specifically ATCAL) in 
the field of aggregated combat modeling will greatly enhance model believability and 
overall usefulness. As computer technologies continue to outpace attrition algorithm 
development, verification, and validation, a continued focus towards efficiency and 
accuracy will ensure the applicability of current and future attrition methodologies. 
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C. INITIATION OF ATTRITION ANALYSIS USING 
JANUS 

The key to initiating a detailed study of attrition and its inherent algorithms 
lies within the generation of high resolution combat simulation data. As such, the 
researcher must carefully consider the process of creating high resolution scenarios. 
The interested researcher should consider the following steps which outline the model 
formulation and data collection process. 



Step 1 


Define the tactical scenario (forces composition, terrain type, 
environmental conditions, and concept of the operation) 
for which a model is desired. 


Step 2 


Construct the Janus simulation scenario and validate its 
representation of the desired interaction between modeled forces. 


Step 3 


Determine an appropriate number (n) of iterations for which to 
execute the Janus simulation. 


Step 4 


Independently execute the Janus simulation n times to develop 
a robust sample of simulation output data from which statistically 
significant killer/victim information can be determined. 


Step 5 


Analyze the consolidated high resolution output data and estimate 
the attrition (kill) rates between opposing weapons systems. 


Step 6 


Incorporate the derived attrition (kill) rates in the analysis of 
attrition algorithms inherent to ATCAL (and other aggregated 
attrition methodologies). 


Step 7 


Utilize an appropriate Lanchester-type deterministic model to 
replay the simulated battle. 



Table XI. Steps for Attrition Analysis and Deterministic Battle Replay Using Janus 
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APPENDIX A. JANUS COMBATANTS 
LISTING AND RUN-TIME BATTLEFIELD 

PICTURES 



This Appendix provides the reader with the detailed description of the combat 
systems making up the REDFOR and BLUEFOR modeled in the AOAC 97-3 Janus 
Simulation described in Chapter II. Additionally, pictures are included to depict the 
battlefield details and force allocations prior to the start of hostilities and at 30 minute 
increments throughout the simulation. 



T-86Main Battle Tank 15 

BMP-2 Armored Personnel Carrier 55 

BRDM-A Armored Personnel Carrier 10 

BRDM-2 Armored Personnel Carrier 7 

GMZ-MI Armored Mine Warfare System 2 
MTKJLI Armored Engineer System 1 

2S12 M 122mm Self-Propelled Artillery 16 

IF TRU Wheeled Truck 26 

DF Tru Wheeled Truck 16 

2S1 155mm Self-Propelled Artillery 30 

BM-21 Rocket Launcher 8 

SNAR-2 Radar 2 

ZOOPAR Radar 2 

2S5 152mm Self-Propelled Artillery 36 

ADA TM Air Defense Artillery Team 5 

ZSU-23 Air Defense Artillery 3 

HIND Attack Helicopter 6 

POL AI Aviation Fuel Depot 2 



Total REDFOR Combatants 242 



Table XII. REDFOR Organization for AOAC 97-3 Janus Simulation 
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M1A1 Main Battle Tank 82 

M1A1P Main Battle Tank w/ Mine Plow 24 

M1A1R Main Battle Tank w/ Mine Roller 8 

M2 Armed Armored Personnel Carrier 113 

BSFV Armed Improved Armored Personnel Carrier 6 
Ml 13 Armored Personnel Carrier 21 

FIST-V Artillery Command Vehicle 4 

CEV Combat Engineer Vehicle 4 

BRIDGE Mobile Bridge Carrier 5 

ADA TE Air Defense Artillery Team 4 

ACE Armored Combat Earthmover 13 

AVLM Armored Vehicle Launched Bridge 9 

AVLM-W Armored Vehicle Launched Bridge 6 

AVLB W Armored Vehicle Launched Bridge 2 

HMV.50 Wheeled Vehicle 47 

HMV/M1 Wheeled Vehicle 10 

HMV/SM Wheeled Vehicle 6 

HMV/CO Wheeled Vehicle 8 

IF TRU Wheeled Truck 25 

DF Tru Wheeled Truck 4 

M106A1 120mm Mortar Vehicles 22 

M577 Armored Personnel Carrier (Command) 6 

GSR/Ml Ground Surveillance Radar 1 

SCT TE Scout Team 34 

SAPPER Combat Engineer Team 90 

SAW TE Machine Gun Team 144 

VOLCANO Mine Dispensing System 4 

M109A3 155mm Self-Propelled Artillery 24 

M109A6 155mm Self-Propelled Artillery 18 

HMMQ36 Artillery Radar System 1 

ADA TE Air Defense Artillery Team 4 

SMK113 Armored Vehicle Smoke Generator 3 

A-10 Fixed Wing Close Air Support 12 



Total BLUEFOR Combatants 764 



Table XIII. BLUEFOR Organization for AOAC 97-3 Janus Simulation 
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JANUS ANALYST 




Figure 1. BLUEFOR Initial Disposition for AOAC 97-3 Janus Simulation 
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Figure 2. REDFOR Initial Disposition for AOAC 97-3 Janus Simulation 
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Figure 3. Battlefield Picture at Start (0 Minutes) for AOAC 97-3 Janus Simulation 
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Figure 4. Battlefield Picture at 30 Minutes for AOAC 97-3 Janus Simulation 
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Figure 5. Battlefield Picture at 60 Minutes for AOAC 97-3 Janus Simulation 
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Figure 6. Battlefield Picture at 90 Minutes for AOAC 97-3 Janus Simulation 
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Figure 7. Battlefield Picture at 120 Minutes for AOAC 97-3 Janus Simulation 
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Figure 8. Battlefield Picture at 150 Minutes for AOAC 97-3 Janus Simulation 
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Figure 9. Battlefield Picture at 180 Minutes for AOAC 97-3 Janus Simulation 
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Figure 10. Battlefield Picture at Completion (211 Minutes) for AOAC 97-3 Janus 
Simulation 
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APPENDIX B. APPLICATION OF THE 
ANDERSON FR METHOD TO A 
HYPOTHETICAL HETEROGENEOUS 
BATTLE REPRESENTATIVE OF THE 
AOAC 97-3 JANUS SIMULATION 

This appendix provides the reader with an example application of the Ander- 
son FR Method with subordinate APP process for attrition rate estimation using 
hypothetical kill rate data representative of the simulation output of the AOAC 97-3 
Janus Simulation described in Chapter II. Specifically, the author has selectively 
chosen to calculate the BLUEFOR and REDFOR Firepower Scores (FPS) and the 
resultant Firepower Index (FPI) using two unique BLUEFOR weapon systems and 
three unique REDFOR weapon systems. The heterogeneous aimed fire battle between 
these opposing weapons systems provides the reader with a numerically obvious closed 
form representation for both forces’ FPS. 

1. FORCES COMPOSITION AND NOTATIONAL TRANS- 
LATION 

Let X represent the BLUEFOR composed of two weapon system types, namely, 

114 Ml Main Battle Tanks and 113 M2 Bradley Fighting Vehicles. Let Y represent 
the REDFOR composed of three weapon system types, namely, 15 T-80 Main Battle 
Tanks, 55 BMP-2 Armored Personnel Carriers, and 17 BRDM Armored Personnel 
Carriers. Recall from Chapter III that the goal of the Anderson FR Method is the 
representation of the heterogeneous force as a single scalar measure of combat power. 
Translating Equation III.l from Chapter III, we have the mathematical expressions 
detailed below by Equations B.l and B.2. 
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FPIx 



(BLUEFOR Firepower Score) 



= ZsfXi 

i=i 

i = integer index where 1 =M1 and 2 =M2 
S * = value (score) of one Ml 

S 2 = value (score) of one M2 

Xi = 114 Mis 

X 2 = 113 M2s 



(B.l) 



3 

FPIy = Y, S J Y i (REDFOR Firepower Score) (B.2) 

i=i 

j = integer index where 1 =T-80, 2 =BMP-2, and 3 =BRDM 
S f = value (score) of one T-80 
S\ = value (score) of one BMP-2 
— value (score) of one BRDM 
Yi = 15 T-80s 
Y 2 = 55 BMP-2s 
y 3 = 17 BRDMs 



2. DETERMINATION OF WEAPON SYSTEM VAL- 
UES 

Recall that the fundamental methods for determining the values of Sf and 
Sj include the analysts’ perception of combat value, historical combat performance, 
measures of weapon firepower and support platform characteristics, and, measures of 
what a weapon can kill. The last method attempts to compute the value (score) for 
a weapon as being proportional to the rate at which it destroys the value of opposing 
enemy systems. This method for determining weapons system scores, called the 
AntiPotential-Potential (APP) method, is the primary focus of this appendix. 
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The formulation of the APP problem relies upon the “rate” at which systems 
destroy opposing systems in battle. These assumed or derived non-negative scalar at- 
trition rates (one for each weapon system firing at every opposing weapon system for 
both forces) account for the highly detailed engagement process described in Chap- 
ter I. While the actual kill rates for each pairwise combination of weapon systems 
resulting from the Janus simulation described in Chapter II could be derived through 
analyzing the killer/victim scoreboard, the author has chosen to use hypothetical kill 
rate values for each pair of interacting weapon systems thus defining kill rate ma- 
trices, Kij and Lji below. This simplification assists the reader with focussing on 
the task at hand, understanding the Anderson FR methodology and the subordinate 
APP process. An example derivation of attrition rates from high resolution output 
is discussed in detail in Appendix C. 
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The matrix Kij represents the rate at which one type i weapon system (1=M1 and 
2=M2) kills type j weapon systems (l=T-80, 2=BMP-2, and 3=BRDM). Analo- 
gously, Lji represents the rate at which one type j weapon system kills type i weapon 
systems. As an illustration, the reader will note that the BLUEFOR (X force) Ml 
Main Battle Tank kills the REDFOR (Y force) BMP-2 Armored Personnel Carrier 
at the rate of Ku = -5. 

The remaining unknowns of Equations B.l and B.2 are the weapon systems 
scores, Sf and Sj. Deriving these individual scores equates to solving Equations B.3 
and B.4 (below) for the unknowns Cx,Cy,S? , and Sj. 



3 



C X S? 


= ZKjsr 

j = 1 


for i = 1, 2 


(B.3) 


CySj 


= ELjiSf 
i = 1 


for j = 1,2,3 


(B.4) 
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Cx = proportionality constant 
Cy = proportionality constant 

So, for example, the value of the BLUEFOR Ml Main Battle Tank (S*) is propor- 
tional to (Cx) the sum total value of every REDFOR system that it kills. The reader 
will note that the first equation above represents two equations (one each for i = 1,2) 
and the second equation represents three equations (one each for j = 1,2,3) for a 
total of five equations. Moreover, this system of five equations has seven unknown 
variables, Cx, Cy , S* , S* , S}' , S% , and S% . A system of linear equations with more 
unknowns than equations defining them is referred to as underdetermined. 

a. Developing the Representative Eigenvalue Prob- 
lem 

To solve the value equations defined by Equations B.3 and B.4, consider the 
following substitutions. 



CySJ = £ LjiSf for j — 1,2,3 
1=1 

<^s] = ^r(£ L iiS?) for j = 1,2,3 
Cy t=i 

Substituting the expression for Sj above into Equation B.3 yields 



CxSf = for . = 1,2 

j-l L/V i = 1 
2 3 

^ CxCySf = ttiKijLjtS?) fori = 1,2, (B.5) 

i=l i=l 

and, applying the same construct to Equation B.4 



3 2 

C x C Y Sj = ZZlLiMjSr) for j = 1,2,3. (B.6) 

j— 1 2 = 1 
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For clarity, Equations B.5 and B.6 can be expressed more compactly us- 
ing vector and matrix notation as depicted by Equations B.7 and B.8. Note that 
£? =1 Lji) = (KL) is a, square matrix of dimension (2x2) and J2j=i E;=i (LjiKij) 

( LK ) is a square matrix of dimension (3x3). 
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.2600 


.3500 


.4000 




.4100 


.5500 




LI< = 








KL = 


.2500 


.3400 


and 


.2000 


.2700 


.3100 














.1400 


.1900 


.2200 



C X C Y S X = (KL)S X 
CxCyS* = (LK)S r 



(B.7) 

(B.S) 



By letting A = CxCy [Ref. 9], the reader will note that Equations B.7 and B.S 
can be written as depicted below and represent a pair of eigenvalue problems which 
will share common eigenvalues. 



A S x = ( KL)S X (B.9) 

AS 1 ' = ( LK )S Y (B.10) 

The reader should note that the matrix products ( KL ) and (LK) are square 
matrices which have only positive entries. Consequently, the Perron-Frobenius The- 
orem states that the largest eigenvalue of a matrix with only positive entries and 
the entries of the eigenvector associated with the largest eigenvalue are positive real 
numbers. See [Ref. 12]. 

b. Solving the Representative Eigenvalue Problem 

The eigenvalue problems represented by Equations B.9 and B.10 can be solved 
by finding the roots of the respective characteristic equations, 



det(KL-XI) = 0 
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<^=^> det( 



) = 0 



,4100 - A .5500 
.2500 .3400 - A 



(B.ll) 



and, 



det( 



det(LK - XI) 

2600 - A .3500 .4000 

.2000 .2700 - A .3100 ) 

.1400 .1900 .2200 -A 



0 

0 



where det is defined as the determinant and / is the identity matrix of the appropriate 



dimension. These characteristic equations are simplified to the following equations 



by evaluating the determinant and collecting like terms. If required, the reader can 
review the methods for evaluating determinants of (2x2) and (3x3) matrices in [Ref. 
12 ]. 



det(KL - XI) = (.4100 - A)(.3400 - A) - (.2500)(.5500) 

= A 2 - .7500A + .0019 

and 

det{LK-XI) = (.2600 — A)[(.2700 — A)(.2200 — A) — (.1900)(.3100)] 

— .3500[(.2000)(.2200 - A) - (.1400)(.3100)] 
+.4000[(.2000)(.1900) - (.1400)(.2700 - A)] 

= -A 3 + .7500A 2 - .0019A 
= -A(A 2 - .7500A + .0019) 

Setting both determinants equal to zero and solving for the values of A which 
fulfill these equalities (using the quadratic formula [Ref. 12]), we find that the first 
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characteristic equation yields roots X\ = .7475 and A 2 = .0025, and, the second 
equation yields A] = .7475, A 2 = .0025, and A 3 = 0. The reader is now able to confirm 
the existence of the largest eigenvalue, A* = .7475, which is both real and positive 
as ensured by the Perron- Fro benius Theroem. Substituting this value, A* = .7475, 
into Equations B.9 and B.10 presents the following eigenvector problems for which 
we seek the solutions. 



(. KL-XI)S x = 0 
(LK-X m I)S Y = 0 

Consider Equation B.12 in isolation. 



(B.12) 

(B.13) 



(KL-\~I)S X = 0 



-.3375 


.5500 




' sf 




0 


.2500 


-.4075 




. S ? . 




0 



As there are infinitely many eigenvectors which solve this equality, the author will 
arbitrarily fix S x = 1 to produce a unique value for S x (and hence, a unique eigen- 
vector). This amounts to normalizing the value of the BLUEFOR Ml Main Battle 
Tank equal to one, which will force all of the weapon systems values to be scaled 
relative to the Ml. The reader should be aware that this scaling method is only one 
from a list of many methods. Further information can be found in [Ref. 9]. Solving 
for the second unknown (the value of the BLUEFOR M2 Bradley Fighting Vehicle), 



the reader can confirm that 5^ = .6134. Therefore, S x = 



1 

6134 



However, the reader should be aware of the fact that this solution vector is 
not unique. Any scalar multiple of this value vector will also satisfy Equation B.12. 
The proportionality constants, Cx and Cy, appearing in Equations B.7 and B.8 will 
simply adjust to absorb the scalar multiple changes to the eigenvector of weapon 
system values - a condition for which the chosen scaling method shortly corrects. 
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Now consider the effectiveness of the BLUEFOR Ml Main Battle Tanks as 



they destroy the “value” of the RED FOR weapon systems. This can be evaluated 

3 

mathematically by the expression ^ KijSj . The IDAGAM scaling method [Ref. 

3 = 1 

3 

9] seeks to scale the values of Sj such that ^ -Ri jSj = Vx*, where A* = .7475 

j = i 

(the largest eigenvalue). Recall that in the development of the eigenvalue problems 
represented by Equations B.9 and B.10 that the author let A = CxCy. Combining 

3 

Equation B.3 with the known information S * = 1 and ^ AT jSj = VxF the reader 

3 = i 

will notice the following result. 



CxSf = jSj 

i=i 

C x = x/A* 

Cx = \[CxCy 
(Cx ) 2 = CxCy 
Cx = Cy = x/F 

3 

Substituting the kill rate information from matrix K into the equation ^2 jSj = 

3 = i 

Vx* results in the equation .45^ + .5 + .5S% = .8646. This particular scaling 
method now calls for the replacement the first equation of the linear system repre- 
sented by Equation B.13 with this new equation. This action results in the following 
matrix equation. 



.4000 

.2000 

.1400 



(LK-X*!)^ = 0 



.5000 


.5000 




sr 




.8646 


-.4775 


.3100 




st 


— 


0 


.1900 


-.5275 




s I 




0 
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Several methods exist for solving this system of three linear equations in three un- 
knowns (many of which are detailed in [Ref. 12] and [Ref. 13]). When solved, the 



resulting values for the REDFOR weapon systems are S Y = 



.8176 
.6309 • 

.4442 



3. CALCULATING THE RESULTANT FIREPOWER 
INDICES 

Armed with the values for each of the BLUEFOR and REDFOR distinct 
weapon systems, the reader can easily calculate the single-valued scalar descriptors 
FPU and FPIy. 



FPU = 



E SfXt 

t= 1 

(1 )(1 14) + (,6134)(113) 



= 183.3142 



FPU = ESjYj 

j-1 

= (.8176)(15) + (.6309)(55) + (,4442)(17) 
= 54.5151 



As described in Chapter III, the Force Ratio (FR^y) can be calculated as FR^y = 

^ J* = l 8 ; 3 f,l 4 : = 3.3626. In a simplified sense, this resultant FR^y would de- 

FPIy 54.5151 

termine the outcome of a battle with respect to parameters set within an aggregated 
model. For a more detailed treatment of the theoretical foundations of the Anderson 
(FR) Method and a critique of its use in force-ratio-based attrition models, see [Ref. 
9]. 

The reader may have gathered that the complexity associated with hetero- 
geneous battles made up of many distinct weapon systems would require the intro- 
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duction of an automated capability to solve the representative eigenvalue problems. 
For example, thoroughly accounting for the 33 BLUEFOR and 18 REDFOR distinct 
weapon systems modeled in the Janus simulation described in Chapter II would result 
in a linear system of 51 equations in 53 unknowns (counting Cx and Cy). The reader 
may also surmise that the kill rate matrices associated with this complex simulation 
may be sparse; in other words, many of the entries of the kill rate matrices could be 
equivalent to zero based upon weapon systems which do not directly interact on the 
battlefield. This fact introduces another consideration to the Anderson FR method 
and the subordinate APP process — the possibility for ill-conditioned matrices (see 
[Ref. 12]). At this point, the reader should have a basic understanding of the An- 
derson (FR) based attrition model and its application to a simplified heterogeneous 
battle. 
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APPENDIX C. APPLICATION OF THE 
COMAN MAXIMUM LIKELIHOOD 
ESTIMATION METHOD TO A 
HOMOGENEOUS BATTLE WITHIN THE 
AOAC 97-3 JANUS SIMULATION 

This appendix provides the reader with an example application of the COMAN 
MLE method for attrition rate estimation using portions of the simulation output 
data from the AOAC 97-3 Janus Simulation described in Chapter II. Specifically, the 
author has chosen to apply this method to the time series of attrition between two 
opposing mechanized infantry weapon systems. The aimed fire battle between the 
BLUEFOR M2 Bradley Fighting Vehicles and the REDFOR BMP-2 Armored Per- 
sonnel Carriers is examined as a homogeneous exchange of direct fires. The primary 
weapon systems engaged in this battle are the TOW anti-tank missile (BLUEFOR) 
and the AT-5 anti-tank missile (REDFOR). As such the derived attrition rate pa- 
rameters will be applicable to the Lanchester Square Law (F|F) for aimed fire war- 
fare. Finally, the author develops a practical application of the differential equations 
(Equations C.l and C.2 detailed below) to validate the derived attrition rate esti- 
mates. This replay is analogous to the replay executed within the ATCAL Phase II 
algorithm. (However, the reader should recall that ATCAL Phase I uses an ad hoc 
method for attrition rate estimation rather than the COMAN MLE method.) The 
data for this analysis was exported from the Janus simulation output file and manipu- 
lated in a commercial spreadsheet software package (not included in this Appendix). 
The data appears in Table XIV on page 104 of this Appendix. However, the in- 
terested reader can obtain copies of the associated spreadsheet files from Professor 
Bard K. Mansager, Department of Mathematics, Naval Postgraduate School. 
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1. NOTATION 



Let X represent the BLUEFOR M2 systems (with TOW missiles) and Y 
the REDFOR BMP-2 systems (with AT-5 missiles). The relative Lanchester Square 
Law equations, modified with the COMAN MLE attrition rate estimates, model the 
change in force size with respect to time. These are Equations C.l and C.2. 



d x 

dt 

dy 

dt 

x{t) 

y(0 

a 

b 



-dy(t) 


with x(0) = xq 


(C.l) 


—bx(t) 


with y(0) = yo 


(C.2) 



size of force X at time t 
size of force V at time t 

derived attrition rate at which BMP-2 firers kill M2 targets 
derived attrition rate at which M2 firers kill BMP-2 targets 



Recall Equations III. 21 and III. 22 on page 44 of Chapter III, as shown below. 



^k— 1 i^k ^k— l) 



Ylk = 1 mk—l(tk tk— l) 

2. DEVELOPMENT OF THE ATTRITION RATE ESTI- 
MATES 

Applying Equations III. 21 and III.22 to the Janus output data found in Ta- 
ble XIV, we have the following: 
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4 

Eit=l n k- l(tk ~ tk-l) 

Total No. of M2 Casualties ( K ) 

E*Li(No. of BMP-2 Surviving at f* — 1) • (Time between kth and k — l sj Casualty) 
17 

469493 seconds 

= 0.00217256 M2 kills per BMP-2 firer per minute; and, 



EaL 1 rn k -i{t k - tk-l) 

Total No. of BMP-2 Casualties ( K ) 

Ef=i(No. of M2 Surviving at tk — 1) • (Time between k t h and k — l st Casualty) 

12 

997883 seconds 

= 0.00072153 BMP-2 kills per M2 firer per minute. 

Intuitively, from the standpoint of classic military operations, the resulting 
attrition rate estimates are logical with respect to the modeled combat scenario. The 
REDFOR is executing a layered defense with fortified positions and established fields 
of fire, while the BLUEFOR is on the attack characterized by weapon system exposure 
and less distinct fields of fire. Therefore, as one would expect, the REDFOR BMP-2 
kills the BLUEFOR M2 at a higher rate than vice versa (all else being equal). This 
provides the analyst with a common sense oriented way of justifying his or her inter- 
mediate results. 

3. DERIVATION OF FORCE LEVEL EQUATIONS FROM 
LANCHESTER’S SQUARE LAW (F|F) EQUATIONS 

From Equations C.l and C.2, the informed reader may see that these ordinary 
differential equations can be solved to obtain expressions for force levels as a function 
of time. Consider the following differential equations: 
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d i - -***> 

t = - Sx(i) 

which, when differentiated with respect to the independent time variable t, yield 



d?x 
dt 2 
d?y_ 
di 9 



a dy 

—a— 

dt 

*dx 

— 6 —. 

dt 



dy dx 

Substituting the values for — and — and using the initial conditions asso- 

dt dt 

dated with Equations C.l and C.2, one has the following two (second order, linear, 
constant coefficient, and homogeneous) ordinary differential equations (ODE). For a 



more detailed explanation of ODEs of this type, consult [Ref. 13 ]. 



d?x 

~dd? 

z(0) 

dx( 0) 

dt 



abx(t) 



x 0 



-ay o 



<P y_ 

dt 2 

y(°) 

<M°) 

dt 



bdy(t) 

y 0 

bx o 



initial condition 1 
initial condition 2 



initial condition 1 
initial condition 2 



(C. 3 ) 



(C. 4 ) 



As such, each second order ODE can be solved using a characteristic equation 
after assuming general exponential solutions of the form x(t) = C\ exp ri< +c 2 exp r2 * 
and y{t ) = C3 exp r3< +c 4 exp r4 *. Considering the second order ODE Equation C .3 and 
its initial conditions in isolation, suppose that x(t) = exp r *, where r is a parameter 
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dx t cPx 2 . 

to be determined. Then it follows that — = rexp and ——r = r exp r . Substituting 

at at z 

these values into Equation C.3 yields 



r 2 exp r * = abexp rt 
(r 2 — ab ) exp r< = 0. 

And, since exp r * > 0 for all values of then (r 2 — ab) must equal zero. Therefore, 
r = ±\ZsS are the roots to the characteristic equation. Note that a and b are both 
positive real numbers implying that ±\Z^S are both positive real numbers. Knowing 
this, the complex and repeated root cases can be eliminated [Ref. 13]. General 
solutions to the second order ODE (Equation C.3) can now be represented by 



x(t) 



C\ exp 



ab)t 



+c 2 exp 



-(Vdb)t 



Examining the initial condition 1, 



z ( 0 ) = xq 



<=> Cl + c 2 = Xo, 



and initial condition 2, 



dx(0) 

dt 



<= 3 - Ci 




-ay 0 



-ay 0 . 



The reader will recognize that the resulting two equations representing the 
initial conditions applied to our general solution form a system of two linear equa- 
tions in two unknowns, namely, Ci and c 2 . Solving this system of equations yields 



1 ^ \ A 1 

Cl = - x 0 - \ -yo and c 2 = - 



Xo + \l -*yo ) . Substituting these values back into 
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the general form of the solutions, Equations C.5 and C.6 result. (Without loss of gen- 
erality, Equation C.6 follows from the work presented for Equation C.5; the interested 
reader is encouraged to verify the result.) 



x(t) = \ \ [ x ° ~ \l jy°) ex P ( v/ab)< + I x 0 + ) exp (Va6)t } (C.5) 



y(t) = ^{[yo-\l TXo ) exp (A/a6)t + ( y 0 + \l tx 0 ) exp (v,aM< 



lb 

V 



1 . 

a 



(C.6) 



As an alternative to the exponential representation of the force size equations 



(dependent on timet), one can make use of the hyperbolic functions cosh = 
exp^ exp 0 

and sinh = to transform Equations C.5 and C.6 into 



exp 0 + exp 0 



x(t) = xo cosh 


t 


— y°\ t sinh 

V o 


-t-i 


y(t) = y 0 cosh 


(^/bbj t 


fb . 

— Xo\j T Sinh 


(yfbbj t 



(0.7) 
(C.8) 

The most significant advantage to this form is the insight provided for the parametric 
relationships that exist within the force level equations themselves. Consider the 
comparison of force size at time t compared to the initial force size for a hypothetical 
force X. This can be accomplished by dividing Equation C.7 by xq, resulting in 



z(t) 

X 0 



= cosh 




ab ) t 



—\ hr sinh 

Xq V b 




ab ) t 



(0.9) 



From Equation C.9, the reader can see that the level of force X appears to be de- 
pendent upon three terms; the initial force ratio ( — ), the relative fire effectiveness 

\x 0 / 



a 



measure y-r j , and what can be called the intensity of combat indicator (^V bbj . In 
a battle to annihilation, these indicator values can be used to generalize simulation 
results [Ref. 9]. 
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Hence, one can now deterministically solve for force sizes, x(t ) and y(t), for 
the desired time t using Equations C.5 and C.6 (or Equations C.7 and C.8) and the 
estimated attrition rate coefficients a and b. This practical extension of Lanchester’s 
Square Law (F|F) affords the analyst a straightforward method for estimating force 
size throughout a battle. The force size results can also be used to verify a and b 
themselves by comparing the numerical values of x(t) and y(t) to the analogous 
high resolution results (force size) at various points throughout the simulation. This 
practical exercise is what remains to be examined in this Appendix. 

4. COMPARISON OF JANUS SIMULATION RESULTS 
TO THE DETERMINISTIC LANCHESTER SQUARE 
LAW RESULTS 

Using the values for a and b derived from the Janus simulation data output, 
one can form the two Lanchester Square Law (F|F) Equations C.10 and C.ll below. 



dx 

dt 


— (0.00217256)y(t) 


with xq = 113 


(c.io) 


dy 

dt 


— (0.00072153)a:(t) 


II 

o 

S* 

‘5 


(C.ll) 



Or, equivalently, using Equations C.5 and C.6 the force size equations can be written 
as 



1 

2 

1 

2 




ab)t 



ab)t 




— (V ab)t 



— (V ab)t 



(C. 12) 



(C.13) 



where 



a = 0.00217256 
b = 0.00072153. 
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Calculating the force size values of x(t) and y(t) for values of t = (0, 15, 30, • • • , 210) 
minutes using Equations C.12 and C.13 above, the reader will note that the deter- 
ministic approximations depicted in Table XV on page 105 of this Appendix are quite 
representative of the Janus simulation results. (Note: The deterministic results for 
x(t ) and y(t) are shown as integer results which were derived by rounding the contin- 
uous solutions to Equations C.12 and C.13). Employing the deterministic Lanchester 
equations to replicate the high resolution simulation results is called “replaying” the 
simulation. This is the role that the ATCAL Phase II algorithm plays as part of the 
aggregated models it supports. In this example, the deterministic Lanchester results 
estimate the casualties achieved during the Janus simulation with fairly high accuracy 
as depicted by Figures 11 and 12 on page 106. Note that the maximum absolute error 
of the approximations is no greater than two weapon systems for both forces. 

The small discrepancies can be attributed to several facts. First, the author ar- 
bitrarily decided to single out two opposing weapon systems from a highly integrated 
heterogeneous combat simulation to demonstrate the mathematics associated with the 
homogeneous form of Lanchester’s Square Law (F|F) and the COMAN MLE attrition 
coefficient estimation process. The reader may have already gathered that the high 
resolution simulation presented in Chapter II infers the need for a heterogeneous-type 
attrition analysis construct; however, tractable formulations do not exist for het- 
erogeneous Lanchester-type equations. So, in presenting the simplistic closed form 
solutions associated with the homogeneous problem, the reader is able to gain an 
appreciation for the mathematics behind attrition computations. Second, examin- 
ing the casualty times listed in Table XIV, one can see that the exchanges of aimed 
fire between the REDFOR BMP-2 (with AT-5) and the BLUEFOR M2 (with TOW) 
were focused around several small battles. Each force shared losses during several 
distinct time periods, i.e., intervals where \tk — tfc-i] are small and follow each other 
closely. The larger values of [tk — depict periods of inactivity during the battle. 
Perhaps distinct COMAN MLE attrition rate estimates a and b could be derived for 
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several time sub-intervals representing their corresponding sub-battles. Lastly, the 
reader should recall that this research only utilizes one stochastic sample of the high 
resolution simulation. A more statistically robust analysis could have been made 
had the attrition rate parameters been estimated after considering the time series 
Markov (Poisson) attrition processes for several independent simulations. 

Even with the perceived shortcomings, the reader will agree that the sampled 
deterministic Lanchester replay of the high resolution simulation is quite good. Ad- 
ditionally, through the application and understanding of this simplified example, the 
reader should have a thorough appreciation for the expanse of this subject matter. 
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Casualty 

Time 


Target 


No. 

Lost 


Firer 


k 


Time Interval No. of M2 
ik L k—\ 


No. of BMP-2 
K-i) 


16.41 


M2 


1 


BMP-2 


1 


16:41 


113 


55 


16.56 


BMP-2 


1 


M2 


2 


0:15 


112 


55 


17.42 


M2 


1 


BMP-2 


3 


0:46 


112 


54 


18.27 


M2 


1 


BMP-2 


4 


0:45 


111 


54 


18.28 


BMP-2 


1 


M2 


5 


0:01 


110 


54 


37.26 


M2 


1 


BMP-2 


6 


18:58 


110 


53 


37.49 


M2 


1 


BMP-2 


7 


0:23 


109 


53 


38.10 


BMP-2 


1 


M2 


8 


0:21 


108 


53 


38.11 


M2 


1 


BMP-2 


9 


0:01 


108 


52 


46.32 


M2 


1 


BMP-2 


10 


8:21 


107 


52 


49.24 


BMP-2 


1 


M2 


11 


2:52 


106 


52 


59.29 


BMP-2 


1 


M2 


12 


10:05 


106 


51 


64.14 


BMP-2 


1 


M2 


13 


4:45 


106 


50 


83.12 


M2 


1 


BMP-2 


14 


18:58 


106 


49 


84.20 


BMP-2 


1 


M2 


15 


1:08 


105 


49 


84.21 


M2 


1 


BMP-2 


16 


0:01 


105 


48 


110.58 


M2 


1 


BMP-2 


17 


2:37 


104 


48 


111.13 


M2 


1 


BMP-2 


18 


0:15 


103 


48 


131.51 


M2 


1 


BMP-2 


19 


20:38 


102 


48 


142.28 


BMP-2 


1 


M2 


20 


10:37 


101 


48 


142.30 


M2 


1 


BMP-2 


21 


0:02 


101 


47 


143.03 


BMP-2 


1 


M2 


22 


0:33 


100 


47 


143.10 


BMP-2 


1 


M2 


23 


0:07 


100 


46 


149.26 


M2 


1 


BMP-2 


24 


6:16 


100 


45 


149.30 


BMP-2 


1 


M2 


25 


0:04 


99 


45 


149.48 


M2 


1 


BMP-2 


26 


0:18 


99 


44 


153.40 


M2 


1 


BMP-2 


27 


3:52 


98 


44 


173.36 


M2 


1 


BMP-2 


28 


19:56 


97 


44 


183.35 


BMP-2 


1 


M2 


29 


9:59 


96 


44 



Table XIV. Janus Simulation Data: Coroner’s Report 

[Ref. 4] 
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t 


Janus x(t) 


Lanchester x(t) 


Janus y(t) 


Lanchester y(t) 


0 


113 


113 


55 


55 


15 


111 


113 


54 


55 


30 


109 


110 


53 


53 


45 


108 


107 


51 


52 


60 


106 


106 


50 


50 


75 


105 


106 


49 


49 


90 


103 


104 


48 


48 


105 


101 


104 


47 


48 


120 


100 


102 


46 


48 


135 


98 


101 


45 


48 


150 


97 


98 


44 


44 


165 


96 


97 


43 


44 


180 


94 


96 


42 


44 


195 


93 


96 


41 


43 


210 


92 


96 


40 


43 



Table XV. Comparison of Janus Attrition Data to Lanchester-Type Replay 
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115 




Figure 11. Graph of Force Size Over Time for BLUEFOR Simulation vs. Lanchester 
Replay Results 




Figure 12. Graph of Force Size Over Time for REDFOR Simulation vs. Lanchester 
Replay Results 
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